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ON THE USE OF DOUBLE FM PULSES FOR 
DETECTION OF TARGETS AND MEASUREMENT OF THEIR PROPERTIES 

by 

Lewis Meier 


ABSTRACT 


This report discusses the use of double FM pulses to measure the sonar path 
length and rate of change of sonar path length (range and range rate in the 
monostatic case) and length and rate of change of length or apparent turning 
rate for a long thin target such as a submarine. Equations are derived for 
the standard deviations of the errors of such measurements. It is shown 
that for coherent processing roof top double FM pulses are preferred while 
for incoherent processing parallel, double FM pulses are preferred. This 
report further discusses the use of two hydrophones with the double FM pulse 
to measure, in addition to the above, bearing and aspect of the target. In 
this case coherent spatial processing is used that is exactly analogous to 
the coherent temporal processing. Again equations for the standard devia¬ 
tions of the errors in these measurements are derived v 


INTRODUCTION 


The major problem in submarine detection is to obtain a large enough signal- 
to-background ratio to reliably distinguish targets from noise. One solution 
is bigger and better transducers, but there is a limit to the amount of sound 
that may be forced into the water and, anyway, at a certain point one becomes 
reverberation limited and beyond this point increased power is of no avail. 
One way out of this difficulty is the use of multiple pulses [1]: if the 
echoes from the pulse are added incoherently, the signal-to-background ratio 
is unchanged but the variances of the background and the signal-plus-back- 
ground are reduced - thus improving detection performance. If the echoes 
from the pulse are added coherently not only are the variances reduced, but 
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the signal-to-background ratio is increased - thus yielding an even greater 
improvement in detection performance. Note, however, that this procedure 
works better in the noise-limited use, since in the reverberation-limited 
case each pulse is in the reverberation of the other pulses. The same 
effect as the use of a multiple pulse can be obtained by the use of multiple 
hydrophones spaced sufficiently far apart that the noise on each is inde¬ 
pendent. In this case, naturally, there are no special problems caused by 
the reverberation limitations. 

Adding echoes coherently is not always a straightforward task. For example, 
adding the ouputs of a linear array of hydrophones will not work unless the 
target is broadside to the array — in general, the hydrophone outputs must 
be appropriately shifted in time to add the echoes coherently. Such diffi¬ 
culties are also an opportunity, because they are usually intimately related 
to some property of the target - such as bearing with respect to the array 
in the example - and thus give the possibility of measuring some target 
property. Even incoherent addition of echoes suffers these difficulties 
and provides these opportunities. For example, if two FM pulses are trans¬ 
mitted, the difference between the times of arrival of the echoes will 
differ from the difference between their transmission times if the path 
length is changing. Since these time-difference changes may exceed the 
length of the matched filtered echoes, they must be taken into account by 
appropriately time shifting the first echo before adding it to the second; 
furthermore, in so doing, the rate of change of path length may be estimated. 

The primary concern in this paper is with pairs of FM pulse echoes obtained 
either by transmitting pairs of pulses or receiving a single pulse on a pair 
of hydrophones. (Note that while in general a pair of hydrophones is talked 
about, in actual practice a pair of arrays - obtained for instance by using 
two arrays or using halves of a single array - would be used.) In this case 
how one echo must be modified to be coherent with the other can be deter¬ 
mined by correlating suitably modified versions of the first echo with the 
second. For example, if the two echoes differ only by a time shift, then 
this time shift may be found by determining the peak of the time correlation 
between the two signals. 

In the case of multiple FM echoes - obtained either by transmitting a pulse 
train or by receiving with an array of hydrophones - a different approach to 
coherent addition must be applied. In essence, this approach is nothing more 
than an extension of classical wide-band beam forming and it will be discussed 
in detail in a sequel document [2]. While the primary focus here is on FM 
pulses, many of the results are equally applicable to any broadband pulse such 
as a short CW pulse. 
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1 BACKGROUND 

A significant effort over a long period of time has been spent at SACLANTCEN 
on the problem of detection and the extraction of measurements. No attempt 
will be made to give an exhaustive list of this work, which can be found in 
[3,4,5]; rather references will be made only to what seems most pertinent 
(some of those documents also contain extensive references ot previous works). 
One fact that distinguishes this work from previous work is that it is based 
on the theory presented by the author in [6]. In the first place, the use of 
this theory clears up some discrepancies in experimental results that will be 
discussed shortly and that were first resolved in [7]. In the second place, 
the use of this theory allows inclusion of true doppler effect rather than 
requiring use of the frequency-shift approximation to the doppler effect, 
which is inadequate for FM pulses used in modern matched-filter sonars. 

Consider the matched-filter output for the echo from an FM pulse. From [6] 
the matched-filtered echo is a modified convolution in time shift x and 
time stretch (doppler) a between the spreading function of the target and 
the ambiguity function of the pulse. The spreading function of a rigid narrow 
target may be represented as a line segment in the r-a plane, as in shown in 
Fig. 1. The various measurement techniques to be described in this paper are 
aimed at eliciting properties of the target-spreading function, which are in 



FIG. 1 MATCHED FILTER OUTPUT FOR FM PULSE AND LINE TARGET 


turn closely related to target properties of interest. In particular, the 
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time shift t q and time stretch a Q of the centre of the spreading function 

of the target are related respectively to the path length l between trans¬ 
mitter and receiver via # the target at an arbitrarily defined zero time and 
to its rate of change l by 


t q = SL/c , 1-a = 


(Eq. 1) 


where c is the speed of sound. (Note that in the monostatic case l = 2r 

and l = 2r , where r is the range and r the range rate.) Furthermore, 

the spread of the spreading function in time shift At is proportional to 
the apparent target length L', while its spread in time stretch Act is 

proportional to the rate of change of apparent target length 

At = 2L'/c , Aa = - 2t'/c . (Eq. 2) 

(L‘ and L 1 will be defined in detail later; in the monostatic case they are 

L cos a and - ctL sina respectively, where L is the target length and 
a its aspect.) The slope of the target spreading function is g , where 


6 = - L'/L , 


(Eq. 3) 


and is closely related to the turning rate of the target. (In the monostatic 
case 6= a tan a .) 

From the results of [6], the effect of making a modified convolution with the 
ambiguity function of the FM pulse may be represented as projecting points in 
the T-cx plane onto the a=l axis by lines with slope -l/(Atj+f Q /k), where 

Aty is the time of the FM transmission relative to the time for which the 

pulse would arrive at the target at t = 0 ; f Q is the centre frequency of 

the FM pulse and k is rate of change frequency in the pulse. This process 
is illustrated in Fig. 1 for the points on the ends of the target spreading 
function. This result makes intuitive sense, since the path length via a 
point with doppler a will change by (l-a)cAt (see [6]) in time At, 
which corresponds to a change in delay of (l-a)At; furthermore, the centre 
frequency of the pulse for this point will be shifted to af Q , which 

corresponds to a time shift of -(a-l)f o /k in the original pulse. 

So far these results do not differ from the results that would be obtained had 
the frequency shift approximation to the doppler effect been used. In fact, 
however, a point on the spreading function does not project to a single point 
on the a = 1 axis. Doppler has two major effects on an FM pulse, only the 
first of which appears in the frequency shift approximation: it shifts the 
centre frequency and it changes the slope in frequency; therefore, when a 
doppler-shifted FM pulse is passed through a filter matched to its original 
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form, energy is lost because the pass band of the filter does not match that 
of its input and the remaining energy is spread in time, both because the 
pass band is finite and because the individual frequency components of the 
input are not shifted by the appropriate phase. These two effects mean that 
both resolution and detection performance are degraded. In this paper the 
results of [6] are used to measure quantitatively the degradation in resolu¬ 
tion and single-echo detection performance and also how measurement extraction 
performance is affected by the resolution obtainable. 


The bandpass effect on resolution dominates for narrow-band FM signals while 
the slope effect dominates for large bandwidth signals. It is shown herein 
that there exists an optimum bandwidth for given centre frequency and pulse 
length that maximizes resolution for a given doppler shift. Kramer [8] has 
investigated the degradation of point-target detection caused by the diminua- 
tion of the peak matched-fi1 ter output resulting from slope mismatch. Herein 
both slope and bandpass mismatch are considered and it is shown that an 
optimum value of bandwidth exists for given centre frequency and pulse 
duration that maximizes peak matched filter output for a given doppler. For 
spread targets, on the other hand, the matched filter output will be averaged 
over an interval approximating the apparent length of the target [1]. In this 
case, it is shown herein that so long as the averaging interval exceeds the 
resolution interval of the pulse, increasing bandwidth improves detection 
performance because the fraction of energy lost is decreased. 


Now consider measurement of target properties by use of a double FM pulse. 
The signal parameters of interest are the centre frequency f , rate of 


change of frequency , duration T^ , and times Aty of the individual 

linear FM pulses, where the index i = 1, 2 is used to indicate the first 
and second pulse in the signal. Derived parameters of importance are the 


bandwidths 


B i ■ 


Ik, |T, 


of the pulses, and the separation AT = At T -At-, 


between the centres of the two pulses, AT’ e (f„ /k, + f /k 9 ) 

o-j i ^2 L 


i 


and 


AT" £ (f 


o 2 /k 2 " f o 1 /k l^ 


In general the 


V s 


and 


T.' s 
i b 


will be 


contrainted by other considerations, leaving the choice of the By 1 s 

available for minimizing doppler dependence and minimizing measurement errors. 
In addition, the choice of AT and the signs of k^ strongly influence 

sensitivity of target-property measurements. (Given f 


T i 


and B 


only the signs of k^ are free in determining AT' and AT" .) The use of 

two identical pulses (for which AT" = o) has the advantage that only one 
matched filter is required and disadvantage that it may be hard to distinguish 
which portion of the matched filter output corresponds to which FM pulse. 


First consider what measurements can be obtained from incoherent processing 
of the echoes from two FM pulses - that is, matched filtering of the echoes, 
magnitude squaring the result, smoothing and measuring the length and delay 
of each echo. Figure 2 shows the matched filter output schematically for 
two FM pulses at times ±AT/2. Results are shown for the cases in which 
the first pulse is a down sweep (k^ < 0) and the second pulse is either a 
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downsweep chosen so that AT" = 0 (f Q /k^ = f Q /k-|) or it is an upsweep 
chosen so that AT'=0 (f„ /k 0 = -f„ A,). A similar diagram can be drawn 

°2 L °i 1 


for an initial upsweep. Several results are immediately apparent from this 
diagram. If AT'=0 the average time shift of the centre of the two echoes 
is equal to t q and the average length in time of the two echoes is equal 

to At. Furthermore, the difference between the centres of the two echoes 
is (1-cx Q ) (AT+AT") while the difference in lengths is -Aa(AT+AT"). Thus 

the sensitivity of these two measurements to (a Q -l) and Act is |aT+AT"|, 

which is much greater when aT‘= 0, and is greatest when the down-up sequence 
is used since the effects of At., and f /k. then coincide; however this 

latter result is of minor significance since f /k- is typically much 
larger than At.. 01 


Plaisant [9,10] proposed measuring the time shift of the centre and the 
length of a matched-fi1tered echo by finding the two time shifts at which 
the smoothed, matched-fi1tered echo power is at one-half its maximum value. 
The centre of the echo is then found by averaging these points and the 
length by differencing them. Furthermore, he developed equations for 
determining the standard deviation of the error in determining these half¬ 
power points by dividing the standard deviation of the smoothed power by the 
mean slope of the smoothed power at these points. Finally, he used these 
results to obtain equations for the standard deviations of the errors in 
measuring range, range rate, and apparent target length. Plaisant's results 
are extended in the present paper in several ways. In the first place, 
measurement of rate of change of apparent target length is considered and, 
in the second place the results are extended to the bistactic case in which 
range and range rate are replaced by path length and its rate of change. 
Furthermore, Plaisant uses the frequency shift approximation and considers 
only the case in which the resolution of the signal is greater than the 
smoothing interval, which in turn is small compared with the apparent target 
length (measured in time). These results are valid provided the resolution 
is determined using the time - stretch model of the doppler effect. In 
addition to this over-resolved case, the present paper also investigates 
the under-resolved case, in which the resolution of the pulses is less than 
the apparent target length but still greater than the smoothing interval 
(i.e. point targets) and also the cases of over- and under-resolved targets 
with no smoothing of the matched filter output. 


Next, consider what measurements can be obtained from coherent processing of 
the echoes from two FM pulses: matched filtering of the echoes, correlating 
the two matched filter outputs in time and frequency, and determining the 
position in time and frequency shift of the peak in the magnitude of the 
correlation. Note that for this correlation to be possible, f and f 

°1 °2 

must have almost the same value and hence may both be taken to be f . In 

analogy with incoherent processing, the time shift of the correlation peak 
is (l-a o )(AT+AT"+6T), where 6T is a correction proportional to the 

frequency shift of the peak. (This correction term, which cannot be 
neglected for large |$|» is found only with the time-stretch model of 
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doppler.) Furthermore, from Fig. 2 the second echo is stretched with 
respect to the first by 

l-e{(AT+AT")/[1+B(AT/2-f o /ki)]}«l-8[AT , 7(l-Bf 0 /k 1 )] , 

since BAT is small compared with 1; therefore the two signals will lose 
correlation for large |B| if AT'=0 (and, hence, AT" is large). This 
result is shown herein by more detailed analysis and has been confirmed by 
unpublished sea trial results. 

One might naively think that this same stretch acting on the carrier 
frequency would cause a corresponding frequency shift of the second echo 
with respect to the first and, hence, that the frequency shift of the 
correlation peak would be 8(AT+AT“)f Q /(1 -B f /k-j). In fact more careful 

analysis carried out herein and confirmed by unpublished sea trial results 
shows that — as first pointed out by the author in [7] — the frequency 
shift is just B[AT+(a Q -l )AT"]f o /(l-B f^k^aB AT f Q . 

Obviously use of AT"=0 will cause little degradation in correlation with 
little loss in sensitivity of the frequency shift of the peak to B ; 
however, it will cause the sensitivity of the time shift of the peak to 
a Q -l to be reduced greatly from |AT+AT"| to |AT|. To retain the sensitivity 

to a Q -l and at the same time not degrade the correlation, more complex 

processing must be used in which the first matched-fiItered echo is not 
only shifted in time and frequency, but also stretched to make it coherent 
with the second matched-fi1tered echo. 


In the above analysis the ordinate of the spreading function is 

9 £ , 

9-r c * where £ P 


a = 1 - ^ * where £ p is the path length via a given point; such a 

spreading function may be termed the temporal spreading function. For 
linear arrays there exists an analogous spreading function in which the 

9 £ 0 

ordinate is a = 1 - -g-t- , where x is the distance along the array from 

an arbitrary origin; such a spreading function may be termed the spatial 
spreading function. (These two spreading functions are just two-dimensional 
projections of a three-dimensional spreading function whose independent 
variables are t, a and a x -) Targets at reasonable ranges and not too near 

end-fire will have a spatial spreading function that may be represented by 
a straight line segment i 

ordinate a replaced by a 


a straight line segment in the x-a plane and hence Fig. 1 applies with the 

X 


a = 1 + cos ip, 
x o 


Aa.. = -2 


Corresponding to a and Aa are 
9L 1 


9x 


(Eq- 4) 


9 £ 


where \p is target bearing relative to the array, since = -cos ip 
(In the monostatic case 9L'/9x = -L sina sim)>/2r, where r 


from the origin of the array to the centre of the target). 

3 x’ 


the spatial target spreading function is c B v , where 


is the range 
The slope of 


= 


9L 

9x 


/ L’ 


(Eq. 5) 
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FIG. 2 MATCHED FILTER OUTPUT FOR DOUBLE FM PULSE AND LINE TARGET 
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and is closely related to the parallax of the target. (In the monostatic 
case 8 x = (simj;/2r)tana where (sin<j;/2r) is the parallax.) 

In analogy with the temporal case, the matched filter output for a single 
echo from a motionless target heard by a single hydrophone located at Ax is 
a modified double convolution of the spatial spreading function with the 
spatial ambiguity function of the FM signal, which again may be represented 
as a projection onto the a x = 1 plane. Figure 1 again applies, with a 

replaced by a x and the slope of the projection equal to -c/x. For a 

moving target three dimensions are required, with the output of the matched 
filter being a double projection onto the a=a x =l line. If two hydrophones 

at x = ±Ax/2 are used, Fig. 2 applies, with the line labelled "initial 
downsweep" applying for the hydrophone at -Ax/2 and that labelled "final- 
upsweep" applying for the hydrophone at Ax/2. In analogy with the temporal 
results, a and B may be measured by correlating the matched filter 
x o x 

outputs from the two hydrophones in time and frequency. The time and 
frequency shift of the peak magnitude of correlation are, analogously, 

(1-a )Ax+(l-a )6T (where 6T is a correction proportional to the 

X _ U X X 

o 

frequency shift of the peak) and 8 X Axf Q /(l-3f Q /k). 

Wiekhorst [li] was the first to suggest measurement of the target aspect 
using the frequency difference between the echoes received at two hydro¬ 
phones - a technique he referred to as space/frequency. In analogy, 

Seynaeve and van Marie [ 12 ] proposed using two pulse separated in time 
to measure apparent turning rate. (In this paper - as has already been 
done - the analogy is used in the reverse direction from temporal to 
spatial processing.) Subsequent to Wiekhorst's original study much effort 
was expended at SACLANTCEN on these techniques, most of which involved use 
of the envelopes of the echoes rather than matched-filter outputs. (The 
envelopes of the echo from a linear FM is an approximation to the frequency 
response of the target expressed as a function of time, while from [6] - 
see below - the fourier transform of the matched-filter output is also a 
measure of the frequency response of the target.) This procedure had the 
advantage of easing computer requirements - an advantage of great importance 
at the time the work was carried out, but of little importance now. It 
suffered the major defect that the frequency shift between the echoes was 
shifted into the time domain hence making it impossible to separate the 
effects of o-l or a -1 from those of B or 8 V - Hug [13] first 

x o 

suggested use of matched filters and correlation in time and frequency to 
obviate this problem in spatial processing. He [14] and Plaisant [9,io] 
made further studies of the use of the double correlation approach in 
temporal processing. 

Plaisant proposed finding the frequency shift of the correlation peak using 
the same technique that he proposed for incoherent processing: namely, to 
find the frequency shifts at which the correlation peak falls to one half 
its peak value and average them. Furthermore, he derived expressions for 
the standard deviation of the error in B by using his previous technique 
of dividing standard deviation by mean slope to determine the error 
standard deviation for measuring the two half-peak frequency shifts. His 
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results are extended in the present paper to measurement of the time shift of 
the correlation peak and determination of its error standard deviation. Also 
treated is the use of the phase of the correlation peak to improve the 
estimate of the time shift between the two echoes. This process is in effect 
an updated version of split-beam processing. Wiekhorst in his original work 
suggested using a double FM pulse as well as two hydrophones with incoherent 
processing between the echoes from the two pulses in addition to coherent 
processing between the two hydrophones for each pulse. In addition to 
investigating this approach and determining the corresponding measurement 
error standard deviations using Plaisant's approach, the present paper also 
treats the use of coherent temporal as well as spatial processing for a two- 
hydrophone, two-pulse combination. 

One final topic must be covered before concluding this background. Figure 1 
shows pictorally how the matched-fiItered echo from a moving, turning, line 
target for a linear FM pulse in a non-dispersive medium is related to the 
spreading function of such a target. In [6] there is a detailed derivation of 
the matched filter output <jr(x,l) that will serve as the basis of the detailed 
analysis contained in the body of this paper. With noise added to the noise- 
free result of [6] and the subscript i added to distinguish the first and 
second pulse, we have 

2ttj f' Q (t-t^ ) 

4>?(t,1) = e 1 1 ^(T-x d _) + n.(x) , (Eq. 6) 

where n(x), is noise, the delays x", and x. depend on target properties, 

a i a i 

and the shifted carrier frequency f^ depends on f Q and target properties. 

The modulation in turn may be expressed in the time domain as the convolu¬ 
tion of a resolution function Yp^ and an impulse response h^ or in the 

frequency domain as the product of a window function ly M and a transfer 
function : 1 

M, . f °° M , , , , M 

V T ) = f y FM. (t " t * a o> h i (x')dx, 

-00 1 


'f-(f) = Tp M (f, a 0 ) H”(f). (Eq. 7) 


The modulation function and the window function, which are fourier transform 
pairs, represent the properties of the pulse, while the impulse response and 
transfer function, which are also fourier transform pairs, represent the 
properties of the target. Note, therefore, that a measure of the frequency 
response of the target is obtained by fourier transforming the matched-filter 
output. 
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The organization of the present paper is as follows: first to be discussed 
are the resolution and single-echo detection capability of a linear FM pulse. 
Next are the incoherent and coherent measurements of target properties using a 
double FM pulse. Third are space/frequency measurements using a pair of 
hydrophones and a double FM pulse. The final topic is how all these results 
can be used in practical algorithms. 


2 RESOLUTION 

This chapter is concerned with the spread of lYp^J 2 in r as a function of 

f Q , T.j, B.j and a Q _. It is convenient to drop the indices and use the 

normalized variables defined in [6] whose definitions are repeated here 
(with a replaced a Q ) for convenience 

x = kT ATT /a Q « kT /fZf^, 


y = (v 1 ) •^o r/a o * {a o _1) ^ 

* -v -V 


(Eq. 8) 


z = kTx . 


It was shown in [6] that in most sonar situations we may use v=0 in expres- 
M 

sions for yp M - a practice that will be followed throughout this document. 

In that case, from [6] it follows that lYp^l 2 is independent of the sign 
of x and y; hence they may be replaced by |x|«B A/f ' and 

,j- - x u 0 

|y| »|o -1| /f Q T. To determine the resolution of |yp M | 2 a measure of 

its spread is required. One such measure is the difference between the 

smallest and largest value of z at which |Yp M 1 2 is some fraction of the 

peak value of |yp M | z over z - terms of x the spread is, from Eq. 8, 

1/B times smaller. Since the spread in terms of x is the quantity of 
interest this will be termed the resolution, while the spread in z divided 
by jx| will be terme d the normalized resolution. Note that the normalized 
resolution is /f Q /T' times the ordinary resolution. 

If the fraction defining the spread in z or x is taken to be 4 /tt 2 «0.4, 
it is well known that for no doppler the resolution (in this case the 

difference between the minimum and maximum x at which 
|Ypn1 2 = [sinc(irBx)] 2 = 4/ir 2 , i.e. x =± 0.5/B) 
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is 1/B and the corresponding normalized resolution is just 1/1xf. For small 

|xy| it was shown in [6], that |Yp M l 2 is stretched by the factor 1/(1-|y/x|); 

hence the normalized resolution is l/(|x|-|y|). On the other hand, it was 

also shown in [6] that, for large |x, y|, |yp M l 2 is rectangular, with width 

2|yx|(1 -1y/xj); hence the normalized resolution is 2jy/x|(|x -|y|). The 
former function decreases with |x| for given |y|, while the latter increases 
with |x| for given |y|; therefore a possible approximation to the normalized 
resolution is the maximum of these two functions. Unfortunately this approx¬ 
imation is not adequate near the best resolution, as may be seen from Fig. 3, 
which plots the actual normalized resolution and this approximation to it as 
a function of |x| for |y| = 0.84. (For f Q =3500 Hz and T=0.5 s this |y| 

corresponds to |ot -1| = |£/c| = 0.02, which thus corresponds to | l | = 60 kn 

or |r| = 30 kn in the monostatic case and, hence, is a reasonable upper bound 
on |ot Q -11 as was pointed out in [6].) 

The problem with the approximation just discussed is the inadequacy of the 
large |xy| approximation. In Appendix A a more suitable approximation to the 
spread for large |xv j is derived that is v alid for smaller |xy|. It is 
l x l /{(|x|-|y|)(1- /l- x| / l(|x|-|y|) 2 |y|])} and it reduces to the above 
approximation when |xy is sufficiently large. The corresponding resolution 
is just this function divided by |x|. A second possible approximation to the 
normalized resolution is the maximum of the previously-given small |xy| 
approximation and this new large jxy| approximation. Fortunately, as may be 
seen from Fig. 3, this approximation is more than adequate - especially 
since the definition of resolution is somewhat fuzzy in any case. 

The minimum spread (corresponding to maximum resolution) occurs when the small 
and large |xy| approximations are equal; this equality occurs when the square 
root in the large |xy| approximation is zero - that is for |x| = x* where 

<**)’-2|y| 

or - since the largest root is obviously the required one - when 


x* = |y| + T/(2|y|) + /I + l/(4|yj 2 ) . (Eq. 9a) 


In terms of the original variables this result implies that the optimum 
bandwidth is 

B *« Iv 1 1 f o + 1/ ( 2 l a o -1 l T ) + ^o /T + ]/ ( 2 l a o ’ 1|T)2 - 

(Eq. 9b) 

Since for any fixed B the spread of lYp^l 2 in x increases with 
|o Q -l|, it follows that the B corresponding to the largest |a -1| of 
interest should be used. 
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M 

For jx| > x* the spread of lYp^l 2 increases slowly; in fact it reaches 

its maximum at |x|=°°, for which the maximum normalized resolution is 2|y|. 

The corresponding ordinary resolution at B=°° is 2|a -1|T. On the other hand, 

M ^ 

for |x|<x* the spread of lYp^l 2 increases rapidly; the normalized resolu¬ 
tion is 2|y| for |x| = |y|+l/2|y|. Correspondingly, the ordinary resolution 
2|« 0 "1|T occurs for B = |cx 0 -l | f Q + 1 /(2[a Q -1 |T). For example, for 

f = 3500 Hz, T = 0.5 s and |ot 0 ~l| = 0.02, B = 217 Hz and the corresponding 

resolution is 6.8 ms — not much worse than the 4.6 ms resolution for zero 
doppler. Furthermore the resolution is 20 ms at B = 120 Hz and B = °° Hz. 


3 SINGLE-ECHO DETECTION 

The detection of point targets is first considered, followed by a consider¬ 
ation of the detection of distributed targets. As in the previous chapter, 
the subscript i will be dropped throughout. 

3.1 Point Targets 

Optimum detection of point targets is based on comparisons of the output power 
of the matched filter with a threshold. In this case the signal power in the 

matched filter output is IYp^1 2 i therefore the key quantity of interest is 

the maximum of |Yp^1 2 * Thus this section is concerned with how this maximum 

varies with doppler and how FM signal parameters should be chosen to minimize 
its falloff with increasing doppler. For this purpose we make the somewhat 
arbitrary assumption that we can tolerate a falloff in the maximum of 

|Yp M | 2 to one half of its maximum at no doppler, which is Iy fm (o.o) 1 2 “ 1 • 

Again it is convenient to use the normalized variables defined in Eq. 8. 

From [6] it is apparent that for fixed a Q the maximum of jYpjyj1 8 over z 
occurs near z=0 (at least for peaks of reasonable magnitude), and, further¬ 
more, |Yp(vfi 2 i s independent of v and an even function of x and y for 
z=0; hence the only remaining variables are |x| and |y|. Let y*(|x|) be 
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M 

the value of |y| that for given |xj yields |yp^(0,a o )| 2 = 0.5. Since a Q is 

very close to 1, by the above criterion, we seek the value of |x| for which 
y*(x) is maximized. Figure 4 is a plot of y* against |x|, from which it is 
apparent that this value of |x| is Jx| = 2.6, for which y*w0.58. (A 
discussion of how Fig. 4 was generated is contained in Appendix B.) In 
terms of the original variables the above result implies that to optimize 
point-target detection B should be close to 

B f = 2.6 /f 0 /T i (Eq. 11) 

and that the corresponding maximum doppler is 

ct f = 1 ± 0.58 / /TJf . (Eq. 12) 

If f = 3500 Hz and T = 0.5 s, B 1 = 218 Hz and |a Q -l| = 0.014, which 
corresponds to |&j = 42 kn or in the monostatic case |r| = 21 kn. 


Distributed Targets 


Generally the resolution of the pulses will be sufficient to over-resolve 
the target. In this case point-target estimation is not optimal and 
a detector very close to the optimum may be obtained by averaging the 
matched-filter output power over an interval before comparing it with a 
threshold (see Heering [l]). The length of the interval, which is not 
critical, should approximate the target length in time. If the target is a 
line target and the scatterers are uncorrelated at separations larger than 

the width of |y FM | 2 in t, Eqs. 6 and 7 imply that the echo power after 

integration is the convolution in x of a function representing the target with 

|yp M | 2 , followed in turn by the convolution of the result with the averaging 

interval. Since the convolution operator is associative, this process is 

equivalent to convolving lYp^l 2 with the averaging integral, followed by 

convolution of the result with the function representing the target. Thus in 

M . 

this case the total energy, rather than the maximum power, of yp M is of 

interest — provided yp^ has sufficient resolution (i.e. is sufficiently 

narrow) that most of its energy is contained within a time interval smaller 
than the averaging interval. 

Now the total energy in yp M is by definition the integral IyfmI 2 over 311 T 

(or equivalently all z). But by Parsival's relationship the integral of 

Iy fm 1 2 over 2 equals the integral of |Tp M | 2 over the normalized frequency 

w = f/B. Under the assumption that the approximation v = 0 is valid, 
from (61 

|Tp M (w,y)| 2 « — (y/x;l)r-L(2w;l-|y/x|); (Eq. 13) 
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therefore 

QO CO 

/| Y p M | 2 dz = /1Tp M 1 2 dw = J ~'-(y/x;l) (1 - |£|) . (Eq. 14) 

-oo -00 


If y*(|x|) is selected to make the integral in Eq. 14 fall to 0.5, as 
was done before with the maximum, it is clear that y* = 0.5|x|. In terms 
of the original variables, this result implies that detection performance 
can be made as large as desired by increasing B - provided that the 
integration interval is large enough. The corresponding maximum doppler is 


Id - 1| = 0.5 B/f . (Eq. 15) 

1 o o 

For example, if f = 3500 Hz, T = 0.5 s and B = 217 Hz, this result 

corresponds to |&|< 93 or, in the monostatic case, |r| < 46.5 kn. For 
this case we have previously seen that for |cx o ~l| < 0.02, the resolution 

is 6.8 ms or better. This is much smaller than a typical integration 
interval of 50 to 100 ms. Furthermore, note that if |a -1| < 0.02, the 

total energy in |Yp M 1 2 over t falls only 32 % (or 1.7 dB) from the 

maximum value at maximum doppler. 


4 INCOHERENT MEASUREMENTS 

Now consider measurement of target properties - this chapter will treat 
incoherent techniques and the next will treat coherent techniques. First, 
the measurement of target properties is reduced to an estimation of the time 
shifts and of the spread of the modulation of the matched filter outputs. 
Second, the use of Plaisant's approach [9,ioj to estimate these quantities 
is studied. Third, the standard deviation and correlations of the errors in 
measuring target properties are related to the like properties of estimation 
errors for the modulation time shifts and spreads. Finally, consideration is 
given to the selection of appropriate signal and signal-processing parameters. 


4.1 Basic measurement principle 

Since most of the target properties of interest are time-varying, a time of 
measurement must be specified. The most convenient time is that at which the 
centre of the signal arrives at the target. The delay of the modulation of 
the matched filter output for the echo from the ith FM pulse is, from [6], 


x d = U + (At T + f /k-)]/c , 

a i c-Z i °i 

«[«. + l (At T> + f 0 .A.j)]/c ; 


(Eq. 16) 

I 

s 
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If x. can be estimated -- say by x. - then £ and £ can be estimated 
a i a i 

respectively, 

A = (T d + ? d ) c/2 - £AT 1 /2 
1 2 

l = (~ d2 " x dx ) c/(AT +AT"), (Eq. 17) 

if, as is assumed from now on, At T = - At T = AT/2. 

I 2 I 1 


As in [9,10] the target is modelled herein as a time segment of uncorrelated 

M 

scatterers of length L ; therefore, from [6], h.(x) is a segment of white 
noise of length 



2L* [1 + -^-(At T + f /k.) - a]/c 
(c-£) 2 'i °i 1 


« 2L 1 [1 + L'(At T + f /k-)] /c . (Eq. 18) 

i i 


The apparent length L' , its rate of change L' and a are defined 
by: 


L' = L[cosa + cos(a-y)]/2, 

• 

L 1 = L[a sina + (a-y) sin(a-y)]/2, 

a = £ L cosa/2L‘ , (Eq. 19) 


where a is the target aspect, y is the bistatic angle, and a and y 
are the rates of change of a and y - all evaluated at t = 0. If 

Atj can be estimated - say by At^ - then L 1 and C' may be estimated 
by 

L* = (AT d + A? d )c/4 - L' AT72 , 

~ 2 1 

L' = (A? d - Ax d ) c/[2(AT + AT")]. (Eq. 20) 

2 1 
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4.2 Plaisant's approach to estimation of modulation delays and spreads 

Now consider estimation of and At^ . Plaisant [9,io] developed 

the following technique: The power in the matched-filter output is smoothed 
by taking a moving average over a suitable interval. For each echo the echo 
plus noise and noise levels in the smoothed, matched-filter output are 
estimated and the beginning and end of the echo estimated by seeing where the 
smoothed echo power crosses a threshold level corresponding to the average of 
these two levels. Finally is estimated as the average of the beginning 

a i 

and end of the echo while Atj is estimated as the difference between the 

beginning and end of the echo. Plaisant does not discuss how the echo plus 
noise and noise levels are to be determined; some approaches will be 
discussed later, for the moment it is assumed that errors in determining these 
levels are negligible. Clearly the average of the two levels is one half the 
echo level in the smoothed matched-filter output; hence the two points for 
which the smoothed echo power cross the threshold levels will be referred to 
as the half-echo-power points. 

Two types of errors may occur in the determination of the half-echo-power 
points: global errors and local errors. A global error occurs when noise 
alone exceeds the threshold level (see 19 ], pp 13 to 16). A basic purpose of 
smoothing the matched filter output power is to decrease the probability of 
global errors. If x threshold on echo signal-to-noise ratio (estimated from 
the echo-plus-noise and noise levels) is set and only those echoes whose 
signal-to-noise ratio exceed the threshold are selected for measuring target 
properties, the probability of global error can be minimized. On the other 
hand, local errors occur because of the variation of the smoothed matched- 
filter power about its expected value. In the remainder of this section the 
standard deviation of such errors will be considered in more detail. In 
general, the standard deviation of the local error in measuring a half-echo- 
power points will be found, as Plaisant [9,10] did, by dividing the standard 
deviation of the smoothed matched filter power at the half-echo-power points 
by the slope of the expected smoothed matched-filter power at these points. 

From Eqs. 6 and 7 the expected value of the smoothed matched-filter power 
P-j (t) is 


Pi(T) = 


s. 

3 1 


(t;T s _/ 2) ® |Yp M .(T,cx 0 ) 


® 

T 


_2 

At 


d i 


l (t-t 


d i ;AT d i / 2 ) a e. 2+ a n 2 


(Eq- 21) 


where T $ is the smoothing interval, a g 2 is the echo energy, o n 2 is 

the noise power density (assumed flat over the frequency band of interest), 
and (x) is convolution in t. Note immediately from this equation that 

T 

Ax. cannot be measured if it is less than T 
d. s. 
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Plaisant implicitly assumes that the spread of |yp^(T,a )| is small compared 

with the minimum of T and Ax d ; i.e. its resolution is greater than 

5 i a i 

they are. For the moment we assume that this is true. The smoothed matched- 
filter output power contains echo and noise. For no doppler, the noise and 
echo were modelled by Plaisant, and are modelled herein, by gaussian white 
noise passed through a filter of bandwidth B i ; hence the standard devia¬ 
tion of the smoothed power at the half-power point can be readily calculated. 

For a non-zero doppler the echo is no longer white noise passed through a 
filter of bandwith B i , but rather through a filter of bandwidth B..-|a 0 -l|f , 

because of the frequency shift caused by the doppler effect. (For this 1 

situation the stretch about the centre frequency may be neglected.) The 
result is that Plaisant's equation for the standard deviation in measurement 
of the time of a half-power point must be modified to 


a 


X 


i 


•k 

T i 


r 


B i T s. 
1 


B iV 


, 1 . 1 
57NT 


2( B i -!v 1 )if 0i ) T * 1 < S/N i) 


(Eq. 22 ) 


where 


T. = min[T , At , ] , 


S i 


S/N, ■ o’ (B, - l V l| f 0 ) y(o n * B,’ AT, T ) 


(Note S/N.J is the smoothed power signal-to-noise ratio.) 

A derivation of this result is contained in Appendix C*. 

M 

Now suppose that the spread of |y FM (x,a )| is still small compared with 
T but not with respect to Ax d - 1 In fact, suppose that Ax rf is 

5>i a i a i 

sufficiently small that (x-x. ; Ax . /2)/Ax. may be considered to be a 

°i a i a i 

6 function in Eq. 21. In this case -following Plaisant's analysis of 
LCW measurement with a non-turning target - the variation in smoothed 


Note that in Plaisant’s work [ 9 , 10 ] . and 0^ are noise powers rather 

than energy and noise power density, respectively, as herein; thus the extra 
factor B.jATj in the denominator of S/N^ . 
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output power is assumed to be caused by noise alone. The result is 


a 


r 


i 


(viv'iy 

B i 2 l Y FM i (0,a o^ 


V 


1 

^7 


i 

57WT 


(Eq. 23) 


As with Eq. 22, a derivation of this result is contained in Appendix C. 

M 

Determining a for At, roughly equal to the spread of |y FM (t ,a )| 2 - 
T i i m i 0 

i.e. the transition between Eqs. 22 and 23 as At. becomes progressively 

u i 

smaller — is much more difficult and will not be done herein. Note that the 

M 

factor before the square root in Eq. 23 is the area under |y fm _(t, a oi )| 2 

divided by its height at t = 0 and hence a measure of its spread; there¬ 
fore this factor is consistent with the first factor in Eq. 22, which is also 
a spread. 

Finally consider the modifications required to Eqs. 22 and 23 as T s becomes 

M 1 

smaller and smaller. For T $ of the order of the spread of |y fm (T,a Q )| 2 , 

the situation is too complicated to be treated herein. 

For T s less than 1/B^, which is certainly less than the spread of 
M ^ 

|y FM (T,a o )| 2 , there is no smoothing of either signal or noise; therefore 

the moving average may be dispensed with - that is T $ set to zero. In 

this case, for At. greater than the spread of [Ypf^ | 2 > the variation in 
a i l 

the power (unsmoothed now) at the half-echo-power points has a chi-squared 
distribution with a mean equal to the noise power plus one half the echo 
power; as a result 



<0,a o > 


°i 


f 


1 1 


4 S/H. (S/N i ) : 


(Eq. 24) 


Note that Eq. 24 retains the same form as Eqs. 22 and 23 except that 
B i T s. is now re Pl aced b y unity. For At^ less than the spread of 


result is 


2 only the noise power contributes to the variation and the 



(Eq. 25) 
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where x* < 0 , 

i 

lYFM i ( T *’ a o >| 2 = \ l Y FH i ( 0,a O ) I' 
S/N i = V lY F M . (O,a ° )|2/0n2 ‘ 


As before, detailed derivations of Eqs. 24 and 25 are given in Appendix C. 

These results are for the noise-limited case. In the reverberation-limited 
case - provided that the reverberation power is constant over an interval 

M 

of length equal to the greatest of T s , Ax^ and the spread of |y'p M | 2 and 

that the spread in doppler of the reverberation is modest (both good 
approximations) — the reverberation can be treated as a zero doppler, non¬ 
turning target of infinite extent. In this case reverberation acts like a 
noise of power density a r 2 /B^ , where a r 2 is the expected reverberation 

power; therefore all the previous equations are valid for the reverberation 
limited case provided that a n 2 is replaced by o r 2 /B^ in the equations for 

S/N i and S/N.! . In the remainder of this paper, although only the noise- 

limited case is referred to, it is to be understood that the reverberation- 

limited case can be treated by this artifice. Obviously B. = a 2 /a 2 

i r n 

is the breakpoint between the noise-and reverberation-1imited cases, with 
the noise-limited case applying for B^ greater than this value. Note 

that from Eq. 22 S/N^ is increasing with B^ for the reverberation- 

limited case and decreasing with B^ for the noise-limited case, hence the 

maximum S/N^ is obtained for B^ = o r 2 /o n 2 ; similar comments apply 

to S/N! . 

As stated earlier, these results are based on the assumption that the 
target is a line segment of independent reflectors whose amplitude have 
identical gaussian distributions. In fact there is evidence that a better 
model of a submarine is obtained if the gaussian distribution is replaced 
by a double-gaussian distribution - i.e. the amplitude is chosen, with 
given probabilities, from one of the two gaussian distributions [ 15 ]. 

This model is consistent with the view what the submarine consists of 
numerous small amplitude reflectors plus a small number of highlights. A 
major difference between the gaussian and double-gaussian distributions is 
that for the same expected echo powers the variance in echo power will be 
greater with the double gaussian. Thus if the gaussian hypothesis is 
replaced by a double-gaussian hypothesis in the analysis leading to Eqs. 22 
and 24, the result will be to increase the magnitude of the first terms 
under the final square root. Similar considerations will apply to the 
coherent processing discussed in the next chapter. 
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Given expressions for the standard deviation in the estimation of the echo- 
half-power points a , the standard deviations a and a. of the 

T ' 7 d i 

estimates x d _ and A~ d found by taking the average of and the differ¬ 
ence between the estimated half-echo-power points, are - assuming no 
correlation between the estimated half-echo-power points - 


a//?. 

i 


o. =a /~2 
AT d i T i 


(Eq. 26) 


The errors in the estimated x^ and Ax,, are uncorrelated since the 

d i d i 

sum and difference of two uncorrelated variables with the same standard 
deviation are uncorrelated. 


4.3 Standard deviations and means of the measurement errors 

Given any technique for estimating x . and Ax d - such as the approach 

a i a i 

of Plaisant just covered - then, the second-order statistical properties of 
the estimates T, £, L' and L* are readily obtainable from o , 

Td i 

a. and p„ . - the correlation between the errors in x H and 

A V T d i AT d i d i 

Ax d - by the use of Eqs. 17 and 20. For simplicity, we assume, as is 
a i 

typically true, that o x = o T = o x and = o At = a. ; then 

d i d 2 d d 1 d 2 d 

the standard deviations of the estimates are, by inspection, respectively 


a = a c /{ 1 + [AT'/ (AT+AT") ] 2 }/2 
£ T d 


= a c /T/|AT+AT"| , 


= a. c /{I + [AT7(AT+AT")] 2 }/8 , 

L' AT d 

s = a C /(/2|AT+AT"I) . 

L*' AT d 


(Eq. 27) 


Furthermore, the correlations between £ and £ and between L and L are 

P_= AT' a j (2 oj . 

££ £ £ 


= AT' or /(2 or ) . 

ll* r t* 


(Eq. 28) 
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Finally if p 

d i d 

and l 1 are zero. 


is zero then all other correlations among £, £, ~L 


The results of the previous paragraph assume that the -rrors in and 

At^ are independent of the errors in and At^ . This is clearly so 


for noise-caused errors; on the other hand, to the extent that the two 
echoes are coherent, it will not be so for errors caused by variations in 
the echo power from its mean value. In the next chapter it will be assumed 
that the two echoes are perfectly coherent - a good approximation when 
parallel pulses (AT" = 0) are used as in coherent processing. When roof¬ 
top pulses (AT 1 = 0) are used, as is typically done for incoherent pro¬ 
cessing, the two echoes may not be coherent for a turning target, but are 
likely to be close to perfectly coherent for non-turning targets. If the two 


echoes are perfectly coherent then a portion of these errors in t 


and 


J 1 


and in At, 


and At, will be identical and hence will not cause 


errors in measurement of £ and L 1 , but will cause doubled errors in 
measurement of £ and L' . In this case Eqs. 22 and 25 should be modified 
by doubling the first term under the final square root when using them to 
determine a p and a. , and eliminating it for and a # 

* L £ I*' 


4.4 Selection of signal and signal-processing parameters 


Consider signal parameters first. To remove the correlations in Eq. 28, it 
is reasonable to set AT' = 0. Let f Q = (f Q + f Q )/2 and T = (T^ + T ^)I2. 


For practical reasons f 


J 1 


r « f and T, 
0 2 0 I 


T , approximately if not 


indeed identically; therefore B-| ~ B 2 ^ B 


, M |2 

Iy F m f « 


hr 


M 

FM, 


lY 


M , 
FM 1 


and AT" a 2 f T/B. Usually AT" is much greater than AT; therefore the 

choice of AT is unimportant. From Eqs. 22 to 27 it is clear with regard to 
measurement of £ and L' that B should be made as large as possible 
provided that the signal-to-background ratio is not degraded too much and 

M 

that the spread of |Yp M I is not made too large. From Eqs. 9b and 10b it 
follows that making T smaller increases the B that minimizes the spread 
of |Tpjvj1 2 at the largest doppler of interest and hence allows a larger B 
to be used without making the spread of jyp M | 2 too large. On the other 
hand,.again from Eqs. 22 to 27 it is clear with regard to measurement of £ 
and L' that B should be made as small as possible provided that the spread 
of |y^| is not made too large, since this makes AT" larger. Increasing 
T also makes AT" larger and at the same time, from Eqs. 9b and 10b, allows a 
smaller B to be used while increasing the spread of lYp^l 2 the least. 
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Thus there are conflicting requirements in signal parameters — small T , 
large B for good £ and L' measurement, and large T , small B for 

good £ and L 1 measurements - that must be traded-off. 


The major signal-processing parameters of importance are the threshold on 
signal-to-noise ratio required to make global errors negligible and 
T.»T « T . With regard to local errors it is clear from Eqs. 22 and 

5 o 1 s ^ 

23 that the lower T is, the better, provided it exceeds the spread of 

M 5 __ 

|-Yp M (T,a o ) | 2 . On the other hand, increasing T $ alleviates the global 

error problem by smoothing the noise, as already mentioned. In fact the 
optimum choice for T s to minimize global errors is Ax^ since, from 

Eq. 22, further increase in T degrades the signal-to-noise ratio. Thus 
selection of T $ is a tradeoff between global and local errors. Since the 

smoothed noise power has a chi-squared distribution if the noise is gaussian, 
it is possible to determine an appropriate signal-to-noise ratio threshold 
that makes the probability of global errors negligible, as will be seen in the 
example below. 


As an example consider a signal with f = f 

°1 °2 


3,500 Hz and T 1 =T 2 =AT=0.5 s 


and suppose the maximum |ot -1} of interest is 0.02. From Ch. 1, the B 

that minimizes the spread of |yp M (x,l±0.02)| 2 is roughly 217 Hz, with a 

corresponding spread of 6.8 ms. If B-j = B^ = 217 Hz and the two FM pulses 

in the signal have opposite slopes, AT' = 0 and |AT"| = 160 (which is 
much larger than AT). For T $ = 20 ms, 2BT s ?=s 8 ; hence the smoothed 

noise power has a chi-squared distribution with eight degrees of freedom. If 
the noise level is normalized to 1, for a signal-to-noise ratio of 4 (6 dB) 
the expected signal-plus-noise level is 5 and the half-signal-power level is 
3. A sample of a chi-squared distribution with eight degrees of freedom 
exceeds three times its mean for 0.229% of the time. If the search for the 
level crossings is carried out over 500 ms — a reasonable interval for a 
target of length not exceeding roughly 100 ms whose rough position is known 
from the detection process — there will be 25 samples and hence a 6% pro¬ 
bability of noise exceeding the level. This probability of global error is 
unacceptably high; however for the slightly higher signal-to-noise ratio 
of 6 (7.8 dB) this probability of global error drops to 0.25%. Of course 
if the noise is not gaussian - as is quite likely at the tails we are 
interested in - these figures will change, but the general form of the 
result is the same. In any case, a signal-to-noise ratio of 6 would seem 
to be a good threshold below which measurement should not be attempted. 


Figure 5 gives a plot of the standard deviation of £ as a function of 
apparent length, assuming incoherent echoes for a signal-to-noise ratio of 6 
and |o -1| equal to 0.00 and 0.02 when the signal described above is used. 

On this figure the ordinary lines represent values computed using Eqs. 22 
and 23, while the dashed lines represent interpolations in the region where 

the spread of lYp^l 2 approximates the apparent target length (in time) and 
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hence neither Eqs. 22 or 23 applies. The standard deviation for £ can be 
found by taking that of T and dividing by AT"/2 = 8 for m/s or by 4 

~ v 

for knots. The measurements L' and L' cannot be made for L' less than 
15 m; for L‘ greater than 15 m their standard deviations are the same as 

the corresponding standard deviations of £ and l . Note that the differ¬ 
ence in standard deviation between no doppler and maximum doppler is small - 

probably less than the accuracy of the results. Figure 5 also gives the 

standard deviations for T and 2£/AT" assuming coherent echoes and no- 
doppler effect and all other conditions the same. As before, for L' > 15m 

the standard deviations for L‘ and L' are the same as the corresponding 

error standard deviations of £ and £ . 


5 COHERENT MEASUREMENTS 

This chapter follows the same general outline as that of the previous. First, 

the measurement of £ and 8 is reduced to an estimation of the position 
of the maximum correlation in time and frequency between the matched filter 
outputs for the echoes from the two FM pulses. Second, the use of Plaisant's 
approach [9,10] for the estimation of this position is studied. Third, the 
standard deviations and correlations of the errors in measuring target 
properties are related to the like properties of estimation errors in the 
position of the correlation peak. Fourth, consideration is given to the 
selection of appropriate signal and signal-processing parameters. Finally, 
incoherent and coherent measurement are compared. 


5.1 Basic measurement principle 

Consider, first, frequency correlation. From [2], 


H^(f) « F{L*C(1-B At-r -a)(f+f )-B(f 0 /k i )(f+^-f )]/(LC)} . (Eq. 29) 

1 i °i u i 2c i 

(As in the previous chapter, c z /(c-£) 2 has been approximated by unity.) 
Since the correlation peak in frequency measures the frequency shift between 
H^(f) and H^f), the form 

H^(f) = (f-Af°) (Eq. 30) 
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is sought. From Eq. 29, Eq. 30 will hold if 


l-0(Atj + -r ~) - 3 

‘l K 1 


Af = 8(AT+AT")f-(f -f ) 

°2 °1 


+ e 


■ J f o 2 +f 0, J Vo, V f o, 
(at4 AT")-V- 8 4 


(Eq. 31) 


MM M 

Because H-| and are being observed through the windows r FM and 

M 1 

Tp M , which are low-pass filters with bandwidths B-j and B^ , f Q -f Q 

must be small compared with B-j and B 2 ; therefore there is no reason not 
to select f = f Q = f Q . In addition, Af should not depend upon f 

because otherwise there will be a stretch as well as a shift between HV 
M 
2 

much larger than AT , AT"«0 . Again there is no reason not to select 

k 0 = k, = k and T, = T 0 = T , in which case r FM = r FM - 

2 


and Ho ; therefore AT+AT" should be close to zero or, since AT" is 
AT , AT"«0 . Again there is no reason 
^ 2 — ^ aim T-j = T 2 = T , in which case T* = ^FM = ^FM * 

B 2 = B-j = B and AT" = 0. Finally a and BAty are negligible with 
respect to 1, but f3f /k is not. With all thesl factors Eq. 31 becomes 


Af° = 3 AT f 0 /(l-ef 0 /k) • (Eq. 32) 

0 MM 

To measure Af° we would like to correlate and Ho , but we have 
D D \ c m 

available only <pand <p 9 . From Eqs. 6 and 7, if segments of and 

M ^ m M 

cj >2 of length Tp centred on x Q are Fourier transformed, and 

can be observed through the window provided by Tp M . For the moment, 

assume that T F is sufficiently long to include the entire echo. The 
• D 

first matched-filter output <j>-T must be advanced by At 0 to align the 

D 

echo contained in it with the echo contained in $2 • Further assume for 
the moment that the noises n^(x) and n 2 (x) are zero and let <t>-j(f) be 
the fourier transform of <^(x+x°-Ax°,1) and be Courier 

transform of ^(t+x 0 ,! ). From Eqs. 6 and 7 
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-2itj 


MO = e 


k(xi -T d ) + f(x d -t 0 +At°)1 

l 0 d 1 d l d l Lm (f-f) 

F FM (f V H l (f V * 


MO = e 


‘ 2lTj [ f o (T d " T d 2 } + f ( T d - T °)l M 

? L 2 2 -^(f-f') H 2( f ' f o) * 


(Eq. 33) 


where from Eq. 26 of [2] 



= T d 2 _T d 2 = ' CA/(c-i)](f 0 /k) • 


It follows from Eqs. 30 and 33 that to align <!>•, and <5>^ in frequency, 

M o ^ 

$1 must be advanced by Af , but from Eq. 17 of [2] and Eq. 8 


J "'-[ f 


(B-lvHf 0 )/2] r”„(f-4f, a 0 ) 


-2irj(a -1 )(2fAf-Af 2 )/k 
= e 0 


J-u [f“Af; (B“| a 0 "lI f 0 )/2]Tp M ( f ) ; 


(Eq. 34) 

therefore, 

-^-[f-^-Af 0 ; (B-|a 0 -l|f 0 )/2] * 2 (f) 

* ‘t’A^ 

e -27Tj{f[(«./c)(AT+2Af 0 /k)-Ax 0 3+ Af°[£/c-(l/c)^ + ——^- - -t°+At 0 ]} 

(B- |a Q -l | f 0 ) / 2] ^'(f-Af 0 ) 

-27rjAf°{[£+£(AT/2-f /k)]/c-x 0 } 
we 

rv [f"f Q ; (B-|a o -l|f 0 )/2]$j(f-Af°) (Eq. 35) 

if 

At 0 = (1/c)(aT + 2Af°/k) . (Eq. 36) 
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Thus if <J?(x,l) is advanced in time by At 0 and in frequency by Af° , 

it will be identical to ^(x,!) over the frequency band in common between 

the two, except for a phase shift. This phase shift may be removed from 

correlation in time and frequency between and by taking the 

magnitude of the correlation. The next section discusses one method of 

estimating the position of the correlation peak and hence At 0 and Af° . 

Given such estimates of At 0 and Af° , 3 and 2. may be estimated using 
Eqs. 32 and 36: 



(Eq. 37) 


In general, of course, At 0 and Af° will contain errors and therefore the 
phase shift A<j>° at Ax° , Af Q will differ greatly from that given in 
Eq. 35. In fact, from Eq. 35 


A4>°«27r ^f 0 (Ax°- Ax°) + Af°{[2+2(AT/2-f Q /k)]/C - x°}j (Eq. 38) 

Clearly if |Ax°-Ax°| < 1/2 f Q so that there are no problems with the 
ambiguity in measuring phase, Eqs. 36 and 38 can be used to obtain an 
improved estimate 2‘ of 2 from an estimate A<t>° of A<J>°: 


c 



+ V T 


-A~f 0 (2/C- 



(Eq. 39) 


In general, the segments of <J>^ and ^ when fourier transformed will be 
approximately centred on the echoes and hence x°«2/c. 


If <(>?(x,l) is advanced in time by Ax° , in frequency by Af° , and 

rotated by A4>° , then added to ^(x,!) and the centre x d and spread 

Ax d of the echo contained in the result estimated (by using, for example, 

the technique of Section 3.2), 2 and L' may be estimated using 

Eqs. 16 and 18: 


l = c f d - £(AT+AT')/2 , 
L' = c Ar d /[2-|'(AT+AT* )1 . 


(Eq. 40) 
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5.2 Plaisant's approach to estimation of the correlation peak 


Consider finding the correlation peak in time and frequency between 4>, 

D 

and $2 as follows: select a time advance At , fourier transform segments 

of length T F containing the echo with the segment for <}>? advanced by At 

r D 1 

with respect to that for tp 2 . Cross-correlate the results, 4>' and 4> , 


in frequency and find the position of the peak by averaging the two 
frequencies at which the magnitude of the cross correlation falls to one 
half of its maximum value. Repeat this process for a series of At 
retaining the position in fequency and the magnitude of the peak for each 
t . Finally find the position of the peak in time by averaging the two 
times at which the series of frequency peaks fall to one half its maximum 

Of course it would have been just as valid to shift <J>? by a series of 

frequency shifts and cross-correlate each in time with g> 2 • In either 

case it is clear that the process of estimating At 0 and Af° is not 
independent; nevertheless it will be assumed the correlation between 
errors in these two estimates may be neglected. 


Now consider the error in estimation of one of the half-peak correlation 
points in frequency. As before, the standard deviation aof this error 

is obtained by dividing the standard deviation of the magnitude of correla¬ 
tion at such a point by the slope of this expected value. Near the 
correlation peak the correlation in frequency between the two matched-filter 
outputs when noise is absent is, except for a rotation, the auto-correlation 
of a segment B-(jot Q -1 | - J01AT)f Q wide in frequency of the transfer 

function. In Appendix D it is shown that the one-half correlation points 
are 0.6/T* distant from the peak and hence will be close to the peak if 
B-(|a -1|-|B|AT)f q T* is large compared with one. In this case the 

magnitudes of the correlation at the one-half correlation points are 
approximately equal and only the variation about these values caused by the 
addition of noise is of interest. Computation of the mean and variance of 
the magnitude of correlation is facilitated by the fact that when the 
standard deviation of a complex random variable is small compared with its 
mean, the mean of the magnitude of the random variable is very close to the 
magnitude of the mean of that variable and the variance of the magnitude of 
the random variable is one-half that of the variance of the random variable 
itself (since it can vary in two directions, but only variation in the 
direction of the mean counts in the magnitude). 

The mean and variance of the correlation in frequency between the two matched 
filter outputs is readily obtained, since the noise in each matched-filter 
outpi.t. is a segment of band-passed white noise and the echo may be approxi- 

M 

mated as such if the spread of lYp^l 2 is small compared with the apparent 
target length. Details are given in Appendix D; the result is 


MV 1 ^ 


2[B-(j V lMB|AT)fjT 


Jjl pTTHTi 

* r B Tp LS/N (S/N) 2 J 


(Eq- 41) 
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where 


T* = min [Ax d , Tp] 

S/N = a|(B-| v l|f 0 )T* /(a 2 B 2 Ax d T p ) 


(Note that S/N is the signal-to-noise ratio in output power of the matched 
filters after smoothing over an interval of Tp) Unlike in the previous 

chapter it is not obvious that errors in the measurements of the two half¬ 
peak-correlation frequencies are uncorrelated; in fact it is shown that 
there is a correlation between the variation in the correlation function at 
these two points but that, to the accuracy of this development, it can be 
neglected. 

The problem of finding the errors in the half-peak-correlation points in time 
is exactly analogous except that time replaces frequency and vice versa. The 
result is 


a 


X 


B- | ot -1 | f 
_ 1 o 1 o 

2[B-(| V l| + |B|AT)f o ] 2 


1 /1 

1*2 , 1 I 

r b Tp 

.S/N (S/N) 2 J 


(Eq. 42) 


Given expressions for the standard deviations in estimation of the half¬ 
peak-correlation frequencies and times, and o^ , the standard 

deviations a.* and a. of the estimates Af° and Ax° are 
Af Ax 

°Af = Op//T » % = 0±//7 . (Eq. 43) 

The errors in the estimates Af° and Ax° are, by assumption, 

uncorrelated. In addition, the phase Acf>° of the double correlation at 

Ax° , Af° can be estimated by taking the arctangent of the ratio of the 
real and imaginary parts of the double correlation function evaluated at 

Ax° , Af° . The error in this phase measurement is shown in Appendix D to 
have a standard deviation of 


c.. 

A<J> 


B-la -1 If. 


B-(|a o -l|+|B|AT)f o 



(Eq. 44) 
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5.3 Measurement error means and standard deviation 

Given any technique for estimating Ax 0 , Af° and A<^° - such as Plaisant's 

approach , just covered — then the standard deviations a 0 , a and o' 

3 £ £ 

of the error in the estimates 8 , 2, and 2 are readily obtainable from 
°Ax ’ °Af and a A4> by use Ec|S ' 37 and The ^ ast mu ^tipiierative 

factors in the equations are for 8 , £ and £' , given there are corrections 
close to unity; therefore errors in these factors may be neglected; hence 


°8 = °Af /(f o AT) ’ °- = a Ax C/AT ; a - = °A<J> c/2Tr f o AT ( Ec l- 45 ) 

An important point to note is the inverse dependence of these standard 
deviations on AT, the time difference between the two pulses of the signal. 

Given any technique for measuring the delay and spread of the combined echoes 
- such as Plaisant's approach, covered in Sect. 3.2 - then the second-order 

statistical properties of the estimates i and L' are readily obtainable 
from o x and o At by use of Eq. 40. Note that if Plaisant's approach is 

a d 

used, the signal-to-noise ratios in Eqs. 22 to 25 should be doubled, because 
of the coherent adding of the two echoes. It would appear from Eq. 40 that 
i contains an error that depends on the error in £ , and similarly that 

L' contains an error that depends on the error in p. However, the error 
in l also causes an error in adding the two echoes; this cancels the 
apparent error in Eq. 40, and similarly for L’ and "g . The error in 
aligning the echoes in time will however cause an error in L' ; hence 


j/ < c °At/ + < 4T /Z - c )/(%/ /2 • 

AT o./(2o l ,) = a AT /l/(a AT ) 2 + (a AT ) 2 . (Eq. 46) 

2 r d 



5.4 Selection of signal and signal-processing parameters 


The two key signal parameters determining performance of coherent measurement 
algorithms are AT and B . Obviously from Eq. 45 the spread between the 
two FM pulses should be made as large as possible as long as coherence 
between the two echoes is maintained. Furthermore, from Eqs. 41 to 44 it is 


33 



SACLANTCEN SM-155 


clear that B should be made as large as possible provided that the signal- 

M 

to-noise ratio is not degraded too much and that the spread of |y | 2 is 

FM 

not made too large. As has already been pointed out in Sect. 3.4, from 
Eq. 9b, use of a smaller T increases the value of B that may be used without 

M 

making the spread of |Yp M 1 2 too large. With regard to the incoherent 

processing used in conjunction with the coherent processing to measure £ 
and l 1 , it is clear from Eq. 46 and the discussion of Sect. 3.4 that 
performance is maximized under the same conditions. Thus for coherent mea¬ 
surements there are no conflicts in choice of signal parameters — small T 

and large B are desirable for £ , £ , L' and B • 

The major signal-processing parameters of importance are Tp and the 

threshold on signal-to-noise ratio required to make global errors negligible. 
With regard to local errors it is clear from Eqs. 41, 42, and 44 that Tp 

should be set equal to At h to minimize the probability of global errors, 

a F 

again because that maximizes the signal-to-noise ratio. Selection of a 
suitable threshold on signal-to-noise ratio required to make global errors 
negligible is not a straightforward analytical problem and is best solved by 
use of simulation and sea trials. 

As an example,, consider a signal with f Q = 3.500 Hz, T = 0.5 and AT = 1 s 
and suppose the maximum |a -1| of interest is 0.02. As already noted, the 

0 M 

B that minimizes the spread of |Yp M 1 2 for |a Q -T| = 0.2 is 217 Hz, with 

a corresponding spread of 6.8 ms. For this signal |AT'| = 160. Take 
Tp = 50 ms as a reasonable average value of At^ . Figure 6a is a plot of 

the standard deviation of B* versus apparent length for a signal-to-noise 
ratio of 6 and the signal just described for |a 0 -l| = 0.0 and 0.02 and 

B = 0.0 and 20 mrad/s (a turning rate a of 20 mrad/s corresponds to a 
complete turn in 5 minutes). For comparison with Fig. 5, note that measur¬ 
ing B by L'/L' gives an error standard deviation of 

a « - 77 } + a. ? 0 /V , 

B l' L' L L' 

which for the case under consideration, coherent echoes, and a target of 

40 m, is 10 mrad/sec. Figure 6b is a plot of the standard deviation of £ 
versus apparent length for the same conditions as in Fig. 6a. For comparison 
with Fig. 5 note that it gives a o, of 0.8 kn for coherent echoes and a 

target of 40 m. 

To use the phase at At 0 , Af° to improve the estimate of £ requires 
that |At°-At°| be less than 1/2 f or that the error in £ be less than 






SACLANTCEN SM-155 


c/2f o AT, which for the example is 0.43 kn, to prevent ambiguities in phase. 

It is clear from Fig. 6b that this will not in general be the case for the 
example. If it were not for problems with ambiguity in phase, in the 

'V 

example the standard deviation of the error in £' would be, from Eqs. 44 
and 45, 0.02 kn for a -1 =6=0; therefore it is indeed unfortunate that 

O 

use may not be made of phase to improve the measurement of £ . 


5.5 Comparison of coherent and incoherent measurements 

Although the examples given in Sects. 4.4 and 5.4 by no means consider 
optimized signal and signal-processing parameter values, the parameter values 
used are not untypical of what might be used in practice. Therefore it is 
possible to derive some general conclusions based on these examples and on 
the analysis that preceeded them. In the first place, incoherent processing 

does relatively better in measuring £ , while coherent processing does 
relatively better in measuring 6 . In the second place, even though £ and 
L' are measured incoherently with coherent processing, there is a roughly 
3 dB gain in signal-to-noise ratio made possible by the coherent addition of 
the echoes from the two pulses prior to this incoherent processing. One 
advantage with the coherent processing is that there is no conflict in choice 

of signal parameters for best measurement of £ and 6 as opposed to £ 
and L' . On the other hand, coherent processing breaks down if the apparent 
target length is too short. 

The major reason that incoherent processing does better than coherent 

processing in measuring £ is the fact that AT" is not zero for incoherent 
processing. A signal with non-zero AT" can be used in coherent processing 
at the cost of increased processing complexity. If such a signal is used, it 
is necessary to correlate the second echo with versions of the first echo that 
are shifted in time and are both shifted in frequency and stretched in 
frequency, where the stretch in frequency is necessary to compensate for the 
stretch in frequency of the envelope of the second echo caused by non-zero 
AT" . The advantage of such a signal is that it can be processed incoher¬ 
ently if L' is too small. The disadvantage of such a signal, in addition to 
the increased processing complexity, is the presence of a conflict in choice 

of signal parameters between those that are best for £ and those that are 
best for the remaining measurements. 


6 SPACE/FREQUENCY MEASUREMENTS 


This chapter follows the same general outline as that of the previous one. 
First, measurement of target bearing and apparent aspect is reduced to 
estimation of the position of the peak in correlation in time and frequency 
between the matched-filter outputs for the echoes from two hydrophones (or 
two separated arrays of hydrophones). Second, use of PIaisant's approach 
[9,io] to the estimation of the position of the correlation peak is studied 


35 







SACLANTCEN SM-155 


Third, target aspect and standard deviations and correlations of the errors 
in bearing measurement are related to the like properties of estimation 
errors in the position of the correlation peak. Fourth, consideration is 
given to selection of appropriate signal and signal-processing parameters. 


6.1 Basic measurement principle 

Consider the geometry given in Fig. 7 and how the matched-filter output for 
an FM pulse depends on the bearing i p with respect to the array, the aspect 
a , and the hydrophone displacement Ax R from an arbitrary reference 

point along the receiver reference line. If the range r R of the target 

from the receiver reference point is large, the change Act in aspect 
caused by Ax H is 

Act = Ax h sin t|yV R . (Eq. 47) 


Now let H^(f,Ax H ) be the H^(f) corresponding to a hydrophone at Ax H . 

The only variable in Eq. 29 significantly changed by the change in aspect 

M 

from a to a + Aa is L 1 ; therefore to find H.(t,Ax u ), L' must be 

3I< 1 M 

replaced by (L'+Aa ) . Furthermore, it is obvious from Eq. 26 of [6] 

and from Eq. 16 that t \ - t . is unchanged but that, in x . , i must 

°i a i a 

be replaced by £+A£ where, from Fig. 7, 

Ai = Ax cos ij> (Eq. 48) 


Now suppose that there are two hydrophones located at ± Ax^ . It is clear 

that differences in L' and i - whether caused by the use of two hydro¬ 
phones separated in space or by two identical pulses separated in time - 
have exactly the same effect; therefore from Sect. 4.1 it is clear that if 
the matched-filter output for a single pulse for the hydrophone at -Ax H is 

correlated with that for the hydrophone at Ax R , the correlation peak will, 

in analogy with Eqs. 32 and 38, have a correlation peak at 

Af ? = 6 x Ax f 0 . f 0 . /k i> * 

Ax° = (Ax cos \p + 2 l Af?/k.j)/c, (Eq. 49) 
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^ |a 0 -l|=0.02, |3|=20 mrad/s 
ao-1| = 0» |B|=20 mrad/s 


FIG. 6a APPARENT TURNING RATE ERROR STANDARD DEVIATION FOR EXAMPLE 



a 0 -l1=0.02, I $1 =20 mrad/s 


S |ot 0 -1 1=0, IeI =20 mrad/s 


FIG. 6b PATH LENGTH RATE ERROR STANDARD DEVIATION FOR EXAMPLE 


ECHO-PLUS-NOISE LEVEL 



FIG. 7 ESTIMATING ECHO-PLUS-NOISE POWER LEVEL AND 
FINDING THE HALF-ECHO-POWER TIMES 
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where Ax - 2Ax^ is the separation between the two hydrophones and - in 
analogy with g - 


3 


x 


_L (ik'Aa + 3L/AX 

1 ^, y9a Ax 3y Ax 


2 sina sin ^ 
cosa + cos(a-y) r R 


(Eq. 50) 


from Eqs. 18 and 47. (Note that from Eq. 
and that the change in y caused by Ax 


3 s - - [,(&: * t ft: |i) 

is the same as that caused in a.) 


These relationships may be used in two measurement schemes based respectively 
on incoherent and coherent time-processing of the matched-filter outputs for 
each hydrophone from the two FM pulses. In the incoherent approach, the two 
matched-filter outputs from the two hydrophones for the echoes Jrom eachjjulse 
are correlated in time and frequency to obtain the estimates At? and Af9 

for i = 1, 2 and then the phases A$^ at these points estimated. Next, 

for each pulse the two matched-filtered echoes from the two hydrophones are 
added coherently by time, frequency, and phase shifting of the matched- 

filtered echoes from the first hydrophone by At° , Af° and A4>° . Now 

the methods of Ch. 3 can be applied to obtain estimates of £,£,L', and 

L' and the estimates cos ~p and 8 X obtained by averaging the two 

estimates of each obtained from 


cos = (c At.. - 2 i Af° /k^)/Ax , 


f s 


?*. = Af o.^- S f o. /k i)^ f o. A *> 


(Eq. 51) 


with 8 = L'/L', which follow from Eq. 48. Alternatively 3 and B x 
may be estimated, following Wiekhorst m] , from 

3 


( Af 2/ f o 2 )-( Af ? /f 0 .) / [(Af°/k 2 )-(Af°/k 1 ) , 


3 . 


Af° Afg 


(l/f 0i k 2)-0/ f o 2 ^) /| (Af°/k 2 )-(Af°/k 1 ) 


A A , 
Uq- 52) 


which also follow from Eq. 48. In any case, if y and r R can be 
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estimated via other means a may be estimated from Eq. 49, cos 41 , and 
~ ? 

0 . For example, in the monostatic case 0 = T tana sin 4 and 0 = a tana ; 
hence a and tana are readily estimated from Eq. 44, cor4 , and t . 

In the coherent approach a parallel double-FM pulse is transmitted and 
temporal and spatial time and frequency shifts are estimated as follows. 

The estimates Ax^ and Af° of the temporal shifts are found by performing, 

for each hydrophone, correlations in time and frequency between the two 
matched-filtered echoes from the two pulses, averaging the two correlation 
functions, and determining the position of the peak and its phase. In a 

similar manner, the estimates Ax° and Af° of the spatial shifts are 

found by performing, for each pulse, correlations in time and frequency 
between the two matched-filtered echoes from the two hydrophones, averaging 
the two correlation functions, and determining the position of the peak and 

its phase. The estimates Ax? and Af? are used to estimate cos 4 and 

e x 

cos 4 = (c Ax° - 2 £ Af°/k)/Ax, 

X x 

0 ; = Af°(l-0 f 0 /k)/f 0 Ax , (Eq. 53) 


which follows from Eq. 48. Finally, by shifting the four matched-filtered 
echoes - two from each hydrophone - appropriately in time, frequency, and 
phase, they may be added coherently, the centre x d and spread Ax d of 

the combined echo estimated, and Eq. 40 used to determine £ and T' . 

In analogy with temporal processing, the estimate of cos 4 obtained via 
Eqs. 51 or 53 may be improved by use of phase information under appropriate 
circumstances. The procedure will be described only for the situation in 
which Eq. 51 is being used as its application; when Eq. 53 is being used the 

procedure differs only trivially. For the improvement to be possible 

|Ax° - Ax°| < 1/2 f is required to prevent phase ambiguities. In analogy 

with Eq. 39, the phase shift A 4 ° in the double correlation at A~° , 

_ _ X X 

< is 

A4° = 27t jf 0 (Ax° - Ax°) + Af° j(£ + Ax cos 4/2 - £ f Q /k)/c - x° 

(Eq” 



Finally, in analogy with Eqs. 39, 44 and 53 may be used to obtain the 
improved estimate 


cos 4 ' 


c 


f Ax 
0 





[' 


(£+ 


£ f /k)/c - 



(Eq. 55) 
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The discussion so far has assumed that only two hydrophones are used. In fact 
in actual practice it is unlikely that only two hydrophones would be used, 
since the signal-to-noise ratio at a single hydrophone would be small. 

Instead, an array of hydrophones would be used. A general discussion of 
space/frequency technique with arrays is beyond the scope of this document; 
however, an obvious way to use an array is to split it into two subarrays, to 
use ordinary beamforming for each subarray, and to treat the beamformed 
outputs as though they were single hydrophone outputs. Obviously performance 
can be improved if the two subarrays are separated in space, as for example in 
a towed array consisting of two active sections separated by an inert section. 

In this regard it is worth noting that the above procedure can be applied in 
the time domain: namely a pulse train of pulses is used, which is divided 
into two sub-pings, each of which is processed by a beamformer that forms 
beams in doppler rather than in bearing. This approach suffers one defect 
not present in the spatial analogue, namely range ambiguity - which, however, 
can be removed by using any of a number of methods discussed by Rihaczek [ 12 ]. 
Alternatively, two identical sub-pings may be used, each of which consists of 
a series of pulses with different centre frequencies. In this case there is 
no range ambiguity, but for detection the echoes in each sub-ping can be added 
only incoherently rather than coherently as is the.case of pulses with 
identical centre frequencies. For measurement of i and 0 , correlation in 
time and frequency are performed for both frequencies and then averaged. The 
net effect is to divided the standard deviation of the measurement error by 
the square root of the number of pulses per sub-ping. From Eqs. 41 and 42 it 
follows that this improvement is almost as good as dividing the signal-to- 
noise ratio by the number of pulses per sub-ping, as would be the case for 
identical centre frequencies. Furthermore, reverberation problems would not be 
as great as when the sub-ping consists of identical pulses. 


6.2 Plaisant's approach to estimation of the correlation peak 

Since the procedure is the same whether the time and frequency shifts have 
temporal or spatial causes, the discussion of Sect. 4.2 need not be repeated. 
It is only necessary to modify Eqs. 41 and 42 to reflect the fact that cos ip 

and 0 x are being measured rather than £ and B and that correlations 

for two pulses exist: 



i 



56) 
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for i = 1, 2 , where T* is defined as before and 


S/N, ( B 1 -| % -l|f 0i ) T*/(o= it d T f ) 


Finally Eqs. 43 and 44 become 



°f./ 




/ /T . 


(Eq. 57) 


c 


A 


A 

V 


i 


B.-lot -1 I f 
i 1 o 1 o ■ 


8 i"( I“o” 11 + 1 6 X I Ax )f 0 . 



for i = 1, 3 . For coherent temporal processing, the subscripts should be 
replaced by x in these equations and the signal-to-noise ratio doubled. 


6.3 Means and standard deviation of the measurement error 

If incoherent temporal processing is used with Eq. 51 to estimate cos ip and 

"V ~ 

P , the errors in the correction terms due to £ and $ may be neglected, 
was done as before, to yield 

X = J Af, /( V x > ’ °cosip. = °At, c/Ax ’ a cos.J,. = <V/ 2 "V X 

A j II 11 1 

(Eq. 58) 


for i = 1, 2. When the two estimates of each are averaged, the standard 
deviations of the error are reduced by 1//2" if f Q & f Q and . 

On the other hand, if coherent temporal processing is used with Eq. 53 to 

estimate cos ^ and , the standard deviations of the error also obey 

(Eq. 58), with the subscript i set equal to x . In this case the signal- 
to-noise ratios used to determine a.- and a. are doubled, which from 

AT x aT x __ 

Eq. 56 again reduces the standard deviations by about 1/ ^(more for low 
signal-to-noise ratios). Note the inverse dependence in Eq. 58 on Ax, 
corresponding to the inverse dependence in Eq. on AT. 

Suppose, instead, that incoherent temporal processing is used with Eq. 52 to 

estimate 6 and B and that f wf » f an d AT' = 0 . These latter 

o 
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assumptions imply that, from the definitions of AT' and AT", 

- f o/ k l = y k 2 = AT "/2 

that, from Eq. 49, 


Af' 


; Af 2 * Af 0 


* e x f o Ax 


and that, from Eqs. 56 and 57, 


°Af 1 « a Af 2 ~ a Af 

Consider first B x , since the analysis is easier. From Eq. 52 and the 
above 


? x »Af°/f 0 Ax (Eq. 59) 


and hence Eq. 58, with the subscripts dropped, applies. On the other hand, 
from Eq. 52 and the above. 


(A~i ■ rf ?) 

2 = AT" 8 Ax f (Eq - 60) 

A U 

and hence 


°6 ~ 2o Af / t AT " 6 x f o Ax ^ - (Eq. 61 ) 


6.4 Selection of signal and signal-processing parameters 

The same factors govern the selection of signal and signal-processing 
parameters for coherent spatial processing as they do for coherent temporal 
processing, with Ax taking the place of AT ; hence the comments of 
Sect. 4.4 apply direct and will not be repeated. Note that if incoherent 
temporal processing is used they will not be repeated. Note that if 
incoherent temporal processing is used there will be a conflict in the 
selection of signal parameters between those that are best for some temporal 
measurement — large T, small B, and those that are best for other temporal 
and spatial measurements - small T, large B. Of course such a conflict will 
not exist if coherent temporal processing if used. 

As an example, consider a pulse with f = 3500 Hz, B = 217 Hz, and 

T = 0.5 s, and two hydrophones placed 100 m apart. If the two hydrophones, 
were 1 m apart then because of the analogy between this chapter and Ch. 4, 
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Fig. 6a for two such pulses 1 s apart in time would apply direct, with 6 

-1 x 
measured in (km) corresponding to B in mrad/s, and Fig. 6b would apply 

with the error standard deviation for cos tp found by taking that for 

l in knots and dividing by two. Since the hydrophones are 100 m apart, from 
Eq. 58 the standard deviations of the error are less by a factor of 100 ; 

therefore B x and its error standard deviation in (km -1 ) are found by taking 

6 and its error standard deviation in mrad/s and dividing by 100 and the 

error standard deviation for cos tp is found by taking that for & in 

knots and dividing by 200. Suppose that ip = 0, r^ = 20 km, and y = 0 ; 

then from Eq. 50 6 = 0.2 km ^ corresponds to tan a = 4 or a = 76°. An 

A 

error standard deviation for 8 X of 0.005 - not unreasonable in the light 

of Fig. 6a — corresponds in this case to that for tan a of 0.1 or at 

a = 0, since the error in a is roughly cos 2 a times the error in tan a , 

an error standard deviation for a of 0.1 radians as 6°. Under the same 
circumstances, an error standard deviation for cos \p of 0.01 - again not 

unreasonable in the light of Fig. 6b - corresponds to that for tp of 

0.05 rad as 3°. Now consider the use of Eq. 52 to measure B . 

It follows from Eqs. 58 and 61, that the error standard deviation for B is 
2/AT" 8 V times that for B Y ; therefore for the pulse and hydrophone pair 

assumed above the error standard deviation for 6 in mrad/s is that of 

Fig. 6a divided by 0.8 B v , with 8 V given in (km -1 ). For a = 45°, ip = 0, 

x x i 

r^ = 20 km, and y = 0, B x = 0.05 (km" ) and, from Fig. 6a, an error standard 

deviation of roughly 12 mrad/s results. As in Ch. 4, for the parameters used 
it is not possible to improve the estimate of cos tp by the use of phase 
information, because of phase ambiguities; however, if that were possible, 
the error would be 0.001 rad ^ 0.06°. 


7 ALGORITHMS 


This chapter looks briefly at how the incoherent and coherent measurement 
techniques discussed in the previous chapters might be used in practice. The 
first step required is estimation of the background power, followed by a 
preliminary single-echo detection. Next, either the incoherent or coherent 
measurement schemes are applied or, in an even more comprehensive scheme, one 
of the space/frequency schemes. Finally, the echoes are combined incoherently 
or coherently, as appropriate, and the final detection decision is made as to 
whether or not a target is present. 
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7.1 Estimation of background power 

One single way of estimating the background level is by smoothing the matched- 
filter output power over a suitable interval with a gap in the centre length 
sufficient to exclude an echo if one occurs there. (If the two FM pulses have 
the same centre frequency, and hence only one matched filter is used, the 
presence of two echoes must be appropriately dealt with by using a long enough 
gap or two gaps.) Alternatively, the noise level may be tracked by a suitable 
recursive filtering algorithm, such as a Kalman-type filter, with data from 
detected targets excluded. 

7.2 Single-echo detection 

Given a running estimate of the noise level using either of the above or any 
other technique, single-echo detection can be carried out in the standard 
manner; that is, by smoothing the matched-filter output power over an 
interval of either T $ or Tp (depending on the type of measurement 

extraction to be used) followed by comparison with a threshold. The purpose 
of the preliminary single-ping detection is twofold: first to reduce the 
amount of data for further processing, and second to ensure that echoes that 
are processed further have adequate signal-to-noise ratio to make probability 
of global measurement errors negligible. For this reason the threshold should 
be set relatively low to ensure a very small probability of measured detection 
at the expense of a relatively high false-alarm rate. In the final detection, 
following combination of echoes, a second threshold will be set to reduce the 
false-alarm rate to acceptable levels. As will be discussed below, it is 
possible to dispense with this preliminary detection in some versions of the 
coherent processing scheme. 

7.3 Incoherent measurement extraction 

The first step in the application of the techniques described in Ch. 3 is an 
estimation of the echo-plus-background-power level and determination of the 
haIf-echo-power times. The accuracy required of the estimate of the echo- 
plus-background levels was not discussed in Ch. 3 and will be discussed only 
briefly here. In the first place the background- and echo-power level 
estimates do not have to be accurate in any absolute sense, but only need to 
be accurate relative to each other. In the second place the acoustic centre 
and acoustic length of a target are only loosely tied to the geometric centre 
and length; therefore a fair error in estimation of the echo-plus-background 
level would seem tolerable. A simple method of estimating the echo-plus- 
background level is as follows: once the smoothed matched-filter power 
exceeds the detection threshold, look ahead by an amount equal to the maximum 
target length plus T $ . If the smoothed matched-filter power at this point 

exceeds the detection threshold the target is too long and the echo should be 
rejected. On the other hand, if this is not true, search back in time to the 
point where the smoothed matched-filter power exceeds the detection threshold. 
Next, average the smoothed-matched-filter output power over the interval 
between the two crossings of the threshold to estimate the echo-plus-signal 
power. The half-echo-power times may now be estimated by averaging the 
background and echo-plus-background levels and searching before and after the 
two detection-level crossings to where the smoothed-matched-filter power drops 
below the half-echo-power level. Finally, the centre and length of the echo 
can be estimated by averaging and differencing the half-echo-power times. 
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This whole procedure is illustrated in Fig. 7. 


If both echoes from a target are detected then measurements of 


a, L 1 and 


L' may be determined from the centres and lengths of the two echoes by using 
the equations of Sect. 3.1. On the other hand, if only one echo is detected 

then (assuming AT’ = 0) estimates of i ± i (AT+ T")/2 and L'±L 1 (AT+AT")/2 
can be obtained by using an obvious modification of the equations of Sect. 3.1. 
If the two FM pulses have different centre frequencies, it will be known to 
which pulse the echoes correspond and hence the appropriate sign that applies 
in these measurements; otherwise this ambiguity is present. Accuracy of the 
measurements may be estimated by use of the estimated echo-plus-background and 
background power levels to estimate the signal-to-noise ratios for the echoes 
and the use of equations of Sects. 3.2 and 3.3 modified in an obvious way to 
reflect the fact that only one echo is present or, if two echoes are present, 
that they have different signal-to-noise ratios. (If the two pulses have the 
same centre frequency it is reasonable to average the two signal-to-noise 
ratios of the echo; otherwise, it is more reasonable to estimate a a_ and a 


At, 


for each echo and average the results.) 


T d 


7.4 Coherent algorithms 

If only one echo is detected in the preliminary detection process described in 
Sect. 7.2, incoherent processing techniques must be used. If two echoes are 
detected, the next step is cross-correlation of the two echoes in time and 
frequency. This task may be accomplished as follows. For the second echo, 
find the maximum in smoothed matched filter output power in the interval Tp 

after the detection time; a segment of the unsmoothed matched filter output 
Tp in length and centred on this maximum is used in the cross-correlation. 

For the first echo, similar segments Tp in length centred AT+At earlier 

are used for the set of Ax over which time correlation is sought. For 
each first-echo segment, the processing is as follows: the first-echo segment 
is multiplied time by time by the second-echo segment and the result fourier 
transformed to yield the cross-correlation in frequency between the two 
segments in frequency. The resulting set of frequency cross-correlations 
is the required time and frequency cross-correlation between the first and 
second echo, which may be processed as described in,Sect. 4.2 to yield the 
time and frequency shifts of the peak, from which l and 6 may be 
estimated by using the equations of Sect. 4.1. 

To estimate l and L' , incoherent processing must be used. To that end, 
segments of suitable length centred on the two echoes are added coherently; 
for this process the segment containing the first echo must be centred on a 

time AT+Ax° earlier than the centre of the segment containing the second 

echo and must be shifted up in frequency by AT° and rotated in phase by 

A 4 >° . It is convenient if the result of the coherent addition is halved so 
that the final result has the same background level as the original echo. 

This result is then smoothed over an interval T , processed as in Sect. 7.3 






SACLANTCEN SM-155 


to find its centre and length, and equations of Sect. 4.1 used to estimate 

Jt and L 1 . Accuracy of the estimates £, 6 , % and L may be obtained 
using the equations of Sects. 3.2, 4.2, and 4.3, taking note of the comments 
made in Sect. 7.3 on the subject of measurement accuracy. 


7.5 Space/frequency algorighms 

For each hydrophone the background level is determined as in Sect. 7.1 
and the two results averaged (unless there is some reason to believe that the 
background levels of the two hydrophones will differ). Next, the single-ping 
detection schemes are applied to the output of each hydrophone. If only one 
echo is detected the incoherent processing of Sect. 7.3 must be used and 
cos ip and $ x cannot be estimated. If only two echoes are detected and 

both are on the same hydrophone, then either incoherent processing as in 
Sect. 7.3 or coherent processing as in Sect. 7.4 is used, as appropriate, 
and again cos t|j and 3 X cannot be estimated. If two echoes are present on 

different hydrophones and it is known that they come from the same pulse then 
incoherent space/frequency processing with suitable and obvious modification 

can be used to obtain a reduced set of measurements (£ and L 1 cannot be 
estimated); otherwise ordinary incoherent processing may be applied to 

each hydrophone and the measurements so obtained averaged (in which case £ , 
L' , cos ip and 6 x cannot be estimated). If three echoes are present and 

it is known which echo is missing, either incoherent or coherent space/ 
frequency processing can be used as appropriate, with, however, suitable and 
obvious modification for the missing echo. If three echoes are present and it 
is not known which echo is missing, the output of each hydrophone must be 
processed individually and the measurements from the two hydrophones 
appropriately averaged (in which case cos <p and B x cannot be estimated). 

The most interesting case is of course when all four echoes — two on each 
hydrophone - are detected. Then.the principles of Ch. 5 can be used to 

obtain estimates of £, £ , L , L' or 3 X • The procedures to be 

carried out should be obvious from the discussions of Ch. 5 and Sects. 7.3 
and 7.4 and will not be belaboured here. 


7.6 Multiple-echo detection 

Multiple-echo detection is in principle no different from single-echo 
detection, except that the echoes are added together either incoherently or 
coherently prior to comparison with a threshold. In the case of incoherent 
processing the matched-filter output power for the first echo suitably shifted 
in time is added to that for the second and then smoothed over an appropriate 
interval prior to the comparison. In the case of coherent processing the 
matched filter output for the first echo suitably shifted in time, frequency, 
and phase is added to that for the second and then magnitude squared and 
smoothed prior to the comparison. If the smoothing interval is Tp , which 

is a reasonable choice, a short cut can be used to obtain the result: add 
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twice the magnitude of the correlation peak squared to the two maxima of the 
smoothed output power used to make the preliminary detection of the two 
echoes. Alternatively, pre-detection can be avoided entirely if the matched- 
filter output power is smoothed over an interval of Tp and at the same 

time a running correlation in time and frequency, such as described in 
Sect. 7.3, is made between the match-filter output at a given time and a 
time AT sooner and added to twice the maximum of the magnitude squared of 
the double correlation for the given time, and comparing it to a threshold. 
Once a detection is made, the techniques of Sect. 7.3 are used to estimate 
target properties. 

Multiple-echo detection for the space/frequency algorithms is similar to that 
just described. If incoherent processing is used, for each pulse in the 
ping, the two echoes from the two hydrophones are added coherently and the 
results for the two pulses are then added incoherently. The spatial analogue 
of the temporal processing just described can be used to perform the coherent 
addition and even to perform preliminary detection for echoes from each pulse. 
If coherent space/frequency processing is used, the four echoes may be added 
coherently by the following spatial/temporal analogue of the temporal 
processing just presented. The matched filter output powers for the two 
hydrophones are added and smoothed over an interval of Tp . 

At the same time a running, temporal, double correlation is made for each 
hydrophone between pulses, the results added, and the maximum magnitude 
squared found. Similarly, running spatial double correlations are made for 
each pulse between hydrophones, summed, and the maximum magnitude squared 
found. Finally, two more running double correlations are made between 
differing pulses on differing hydrophones and the maximum magnitudes found. 

The smoothed power in the coherent sum of the four echoes at a given time is 
the sum of the smoothed, summed, matched-filter output powers at that time, 
the same at a delay of AT, and twice the sum of all the double correlation 
maxima. 


CONCLUSIONS 


With the use of a double FM pulse it is possible to measure £, £, L and L 
or 6 for a long, thin target such as a submarine. These quantities may be 
measured using incoherent processing of the echoes, in which case the best 
signal to transmit is the roof top FM where f /k, » -f /k ? ; further- 

o-| i c. 


more, processing is eased if f is taken different from f so that 

°1 °2 

the two echoes may be clearly identified. In general, to measure £ and L 
well a relatively short pulse of large bandwidth should be used; while to 


measure £ and L well a relatively long pulse of small bandwidth is 


preferable. The quantities £ and L can be measured using coherent 
processing of the echoes, in which case the best signal to transmit 
consists of two identical FM pulses of relatively short time and long . 
bandwidth. Incoherent processing does relatively better in measuring £ , 
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while coherent processing does better in measuring 6 . Use of phase 

information in the coherent processing can significantly improve l 
measurement if phase ambiguities can be avoided. With the use of a double FM 
pulse and two hydrophones it is possible to measure, in addition, cos \p and 
6 X . The temporal processing may be either incoherent or coherent, but 

spatial processing must be coherent. The best signal to use is determined by 
the type of temporal processing to be used. Again, use of phase information 
presents the possibility of significantly improving cos ip measurement if 
phase ambiguities can be avoided. By incoherent and coherent combination of 
the echoes received on either one or two hydrophones, detection performance 
may be significantly improved. 

The analysis of this report assumes that an ordinary, linear FM pulse and only 
a filter matched to it are used. As a result, 4 there is a limit on the 
resolution attainable for targets with large |£| . This limitation may be 
avoided in one of two ways: use of the doppler-invariant FM signal or use of 
a bank of matched filters. With either approach, resolution may be improved 
at will be using a wider bandwidth signal. These procedures might be of use 
if coherent processing is to be employed, as performance may then be limited 
by bandwidth. Most of the analysis contained herein applies if either of the 
resolution-improving schemes is employed, provided that the correct expression 
for resolution is used; however, some of the doppler correction terms may be 
modified or eliminated. 

Much work remains to be done in this area. In the first place the equations 
contained herein for the standard deviations of the measurement error should 
be verified by simulation and sea trials. (Some verification of the measure¬ 
ment of Z with sea-trial data has been carried out by Plaisant [ 10 ], who 
indicates that his equations and the extensions contained herein may be 
optimistic.) A second area of major interest is the effect of time and 
frequency spreading on measurement errors. Finally, the question of what the 
measurement procedures will yield for targets that are notlong, thin, and 
rigid is of interest since false targets may in fact not be of that form. The 
theory contained in [6] provides the basis for analyzing the effects of the 
medium and the target on the measurement processes. 

One type of dispersive medium that can be analyzed with little further work is 
a medium that spreads the signal in time only. In this case the echo is 
simply what would be obtained without the spreading, convolved with the time¬ 
spreading of the medium. For example, suppose that the medium contains 
multi paths, then if coherent temporal processing is employed, the double 
correlation must be convolved with the autocorrelation function of the 
multipath structure along the time axis. If the multipath delays are large 
compared with the width of the double correlation peak in time, the techniques 
described herein may be used to determine the autocorrelation function of the 
multipath structure. If the multipath structure could be inferred from its 
autocorrelation, then incoherent processing to measure l and L and 
detection performance could be greatly improved by convolution of the matched- 
filter outputs with the multipath structure. Unfortunately, a given auto¬ 
correlation function will not correspond to a unique multipath structure. One 
possibility is to convolve the matched-filter output with all multipath 
structures consistent with the measured autocorrelation of the multipath 
structure to see which gives maximum output. These processes seem worthy of 
further study. 
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APPENDIX A 

SPREAD OF DOPPLER AMBIGUITY FUNCTION IN TIME 


The purpose of this appendix is to find a reasonable approximation to the 
spread of |yp M | 2 in z that is valid for large |xy| . The fourier 

transformation of |Yr M | 2 with respect to z was given in [6] of the 
Main Text as 

i^(y/x;l) n (w, l-|y/x| );(l-|y/x|-|w|) sinc[2Trxy(l-|y/x|-|w| )w] . 

(Eq. A.l) 


If the term |w| could be dropped from the sine function and the factor 
proceding it, inverse transformation would be easy, and in fact dropping 
|w] leads to the inadequate approximation for the spread of |yp M | 2 for 

large |xy| given in [6]. When |w| is dropped, |Yp M 1 2 is a rectangular 

function whose width is the reciprocal of the value of w that yields it 
for the argument of the sine function. To obtain a better approximation it 

is assumed that the spread of |Ypjvj1 2 is the reciprocal of the value of 

w that yields tt for the argument of the sine function - this time 
including one half of the |w| term (inclusion of the full |w| term leads 
to an unacceptable approximation). The value of w sought must obey 
(sine w is positive). 

\ w 2 -(1 -1y/xj) w + 1 /(21xy[) (Eq. A.2) 

therefore it is 


w = 1 -j y/x| - / (l-|y/x|) z - 1/1xy| 


= (|x|-|y|)(l-/l-{x|/C(|x|-|y|) z |y|J/|x| (Eq- A. 3) 


since the smallest root is obviously the one sought. Inversion of this 
result yields the approximation given in the body of the text. 
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APPENDIX B 

DETERMINATION OF THE DOPPLER SHIFT CAUSING THE PEAK MATCHED-FILTER 
OUTPUT POWER TO FALL BY HALF 


Figure 4 of the main text was generated as follows: From Eq. 14 of [6] of 
the main text it follows that we are seeking the value of |y| for which 


C(b) 2 + S(b) 2 = 0.5|xy| 

b = /Ixyf (l-|y/x|) (Eq. B.l) 


since for z = 0, a = -b and C( ) and S( ) are odd functions. By 
selecting values of b and solving these equations, y*(|x|) is readily 
generated. 
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APPENDIX C 


STANDARD DEVIATION OF THE HALF-ECHO-POWER POINTS 


This appendix is concerned with Plaisant's approach to determining the time 
shift and spread of the target echo for an FM pulse by determining the 
times at which smoothed matched-filter output power crosses a level half 
way between the full-signal and no-signal levels. Four cases are to be con¬ 
sidered, depending on the relationship between the apparent target length 

(measured in time), the averaging interval, and the spread of |Yp M J 2 - * n 

all four cases the standard deviation of the error in measuring this time 
is found by dividing the standard deviation of the smoothed-matched-fi1 ter 
output power at this time by the slope of the expected smoothed-matched- 
fi Iter output power at this time. 


C.l Case 1 Spread of !Yp M 1 2 less than both apparent target length 
and averaging interval 

In this case the convolution with |Yp M 1 2 in Eq. 21 may be removed and 

replaced by its integral over t , which is readily found: from Eq. 14 

the integral of |Yp M 1 2 over z is (B^-1« 0 ~11f_ )/B^ ; the integral 

i °i 

over t is just times smaller, or (B^-|a Q -1|f Q )/B? . The signal- 

to-noise ratio is found by taking the ratio of the maximum of Eq. 21 to 
its minimum and is, since the maximum of Eq. 21 occurs when the target is 
centred in the smoothing interval. 


f (B 1 - |c o -i | f 0< )/Bj ]/r r >/^ 


(Eq. C.l) 


where T**is the larger of T and At. . This result is identical to 
1 s i °i 

that given in the main body after Eq. 22 and is valid also aS-Jong as at 
least either T $ or Ax d _ is larger than the spread of lyj^ | 2 - i.e. 

also for the next two cases. Since the expected smoothed echo power is. 
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from Eq. 21, equal to 


(Eq. C.2) 


times the convolution between two j~»- functions, its derivative at the - 
half-echo power point is just A. 

Calculation of the standard deviation of the half-signal-power point is 
harder. The smoothed power at this point consists of 1/T times the 

** "j 

sum of two uncorrelated signals (see [9] of the main text): the integral 
over T*/2 of echo-plus-noise and the integral of T s -T*/2 of noise 

alone. Because the noise is a band-limited, gaussian, white noise with 
power a 2 and bandwidth B i , the latter integral is proportional to a 

chi-squared variable with 2 B.(T s -T*/2) degrees of freedom; hence 

its variance is cr 2 (T s -T|/2)/B i . For no doppler the first inte- i 

gral is also readily evaluated because the echo-plus-noise is again 
band-limited, gaussian, white noise with bandwidth B. , but with power 

o 2 +c 2 /B-Ax , ; the result is a variance of (a 2 +a 2 /B. At. ) 2 (T*/2)/B.. 

n e i 1 a i i i 

Adding these two, taking the square root, and dividing by T , yields a 

i 

standard deviation of smoothed power at the half-signal power point of 




A » 


T «, 4 \ 


i 


i 


2Bi 


^i 


C(°e./Bi 


^ d .) 2 


T ? 


2(a 2 / Bi AT d . 

1 l 


°n T ? + 


K 


T s .] 


(Eq. C.3) 


which is Plaisant's result. 


Unfortunately, for non-zero doppler shift, the echo no longer has band¬ 
width B. and the first integral is more difficult to compute. What is 

required is the variance of 


T*/2 T*/2 T*/2 T*/2 

/ |(t) + X 2 (x)| 2 dT = f 0 |x^( t)| 2 dx + f 0 2Re[x*(x)x 2 (x)]dx + f Q |x z (x)| 2 dx 


(Eq. C.4) 


I 
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where 


V 1 


indicates "real part of" and x 


and x are independent, 

2 


band-limited, gaussian noises of powers o and 


■V B rl 


■'l f o > B i 


Ax 


d i 


B i 


1“ -l|f„ 

1 o 1 o- 


respectively. Because 
the three integrals on the right 


respectively and bandwidths B. and 

of the independence of x and x , 

1 2 

hand of Eq. C.4 are uncorrelated; hence the variance of the integral on 
the left is the sum of the variances of the three on the right. The 
variances of the first and third integrals are readily found to be 

°n T i/( 2B i) and 0 e.( B r 2a o -1 l f o.) T */( 2B i AT d.) » respectively. Since 

x (x) and x (x) are independent, the expected value of the second 

integral is zero and its variance equals the expected value of its 
magnitude squared i.e. 


T*/2 

T i/2 



/ 

/ E{[x*(x)x (x)+x (x)x*(x)]*[x*(o)x (o)+x (o)x*(a)]}dxdo 1 

- T i/2 

- T t/2 

12 12 

12 12 

T*/2 

Tt/2 



=2 f 

/ 

Re E[x (x)x*(a)] ECx* 

k (x)x (a)]dxdo 

-Tt/2 

-Tt/2 

i i 

2 



TT/2 

Tt/2 

=2[V o 2 
n e.j 

(B - - | a 
v 1 1 0 

-l|f 0 ) /B i AT d ] f 
°i 1 a i -Tt/2 

/ sinc[ir(B.-|a -l|f )(x-a)] 

-Tt/2 i 




sincfnB^(x-o)]dxda 

s 2[o^ o 2 
n e^ 

< B cl% 

■ 1 l f 0 . )T T /B f 4T dT r 
1 1-00 

sinc[ir(B 1 .-|a Q -l |f Q )x] sinc(irB^x)dx 

= 2 °n 

n e. 

(®rl“o 

-l|f 0 .) T?/81 4T d . , 

[Eq. C.5] 
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which, when added to the two other portions of the echo-plus-noise variance 
and the noise-only variance, yields for the standard deviation at the half¬ 
signal power point: 



B.-|a -1|f ' 
1 1 0 1 0 • 


(°e. /B i Ax d.> 
i i 


B i 


|T i 


/VK-nv 


+ 2 °n(° ei /B i 


T* + 2 °n 
i n s^ 


J (Eq. C.6) 


Division of this result by A and use of the definition of S/N yields 
(Eq. 22). 


C.2 Case 2 Spread of |Yp M | 2 more than apparent target length but 
less than averaging interval 

In this case the factor in the double convolution (Eq. 21) corresponding to 
the target can be dropped; therefore the slope at the half signal-power- 
point is 



T 


s 


i 





(Eq. C.7) 


Since the variance is due to noise only, the standard deviation of smoothed 
power is readily determined to be 



(Eq. C.8) 


Division of the latter by the former and use of the definition of signal- 
to-noise ratio yields (Eq. 23). 
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M 

C.3 Case 3 Spread of |Yp M l 2 less than apparent target length and 
no smoothing 

In this case the factor in the double convolution (Eq. 21) corresponding to 
the moving average smoothing can be dropped; therefore the slope at the 
half-signal-power point is 





(Eq. C.9) 


Since T is assumed to be less than 1/B- , the variation in smoothed 
s i 

power about its mean is proportional to a chi-squared variable with two 
degrees of freedom. At the half-signal-power point the expected echo power 
is, by definition, one half of the maximum expected echo power; hence the 
standard deviation of the smoothed power is 


(1/2 S/N. + 1) a 2 


(Eq. C.10) 


Division of this result by the slope and use of the definition of signal- 
to-noise ratio yields (Eq. 24). 


M 

C.4 Case 4 Spread of |Yp M 1 2 more than apparent target length and 
no smoothing 

In this case both the factor in the double convolution (Eq. 21) corres¬ 
ponding to the moving average and that corresponding to the target may be 
dropped. The peak signal-to-noise ratio is then S/N. , as given following 

(Eq. 25), and the half echo-power points occur when jyp M (t, a Q )| 2 equals 
one-half. Let t* be the negative value of t for which this occurs; the 

1 j^| 

corresponding slope is just a 2 times the derivative of |yc M (t, a )| 2 

e i r "i 0 

at t| . Furthermore, the variation of the power due to noise alone is 

just o 2 . Use of these results yields (Eq. 25). 
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APPENDIX D 

STANDARD DEVIATION OF THE 50% CORRELATION POINTS 


This appendix is concerned with Plaisant's approach to determining the 
position of the correlation peak in time and frequency between the two 
echoes from a double FM pulse by determining the times and frequencies at 
which the magnitude of correlation is half of the maximum. In analogy 
with Appendix C, standard deviations of the half-correlation-peak points 
are determined by dividing the standard deviation in the magnitude of 
correlation by the slope of the expected magnitude of correlation at the 
half-peak point on the expected magnitude of correlation. Correlation 
between the time and frequency Doints is ignored. In addition, measurement 
of phase at the peak is treated. 


D.l Correlation in frequency 

From the discussion of Sect. 4.1 of the main text it is clear that, for 

Ax - At 0 , the correlation in frequency, if normalized and shifted in 
frequency so the peak is at the origin, is the correlation Y(f) between 
X(f) = y’Wr)] and X ' ( f ) = 9T[x'(t)], where 

x(t) = [Y(T)®h(t) + n(x)]j l(t; T p /2)/ /T^ , 

x'(t) = [y 1 (x) (x) h(x) + n 1 (t)]j—l(t ; T f /2)//T^ , ( Eq Q ^ 

h(x) is a segment of gaussian white noise of length At^ centred on the 
origin and of power At. ; n(T), and n'(x) are uncorrelated, gaussian, 

G Q |y| 

white noises of power a* ; y'(t) is Yp^(T,a 0 ) ; and y(t) is the inverse 

fourier transform of Tp M (f,a o ) with its frequency limits shifted up by Af 0 . 

Thus |E[Y(f)]| a E[|Y(f)|] is sought to determine the frequency and slooe 
of the half-correlation-peak points, var [Y(f)]/2 « var |Y(f)| is souoht 
to find the variance at the half-correlation-peak, and 

cov [Y*(-f) Y(f)]/2 * cov(|Y(-f)| |Y(f)|) 

to find their correlations. 
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The inverse fourier transform of Y(f) is 

y(x) = x*(t) x' ( t ) ; (Eq. D.2) 

therefore, since n(x) and n'(T) are independent 


E[y(x)] = E{[y(x) (x) h(T)]*[Y'(i) (x)h(x)]}- r -T-(x,T F /2)/Tp 

T T 

OO CO 

= j-u(t;T f /2)/T f / / y*(x-c)y‘ (x-cr 1 )h*(o)h(o') dodo 1 


J t -(x‘>Tr/2) AT d/2 


AT d T F 


/ y*(t-o) y'(x-o) da 

_Ax d/2 


— e 


a 2 j— l(t;T*/2) 


"AtTTZ- 

d V 


f y*(a) y'(o) da 


a 2 o -r^(j,J*/2) - 


At , Tj- 

d F -«> 


/ r*(f) r'(f) df 


°g J— ^ (x * T / 2) 


B Aij T f 

^(B-lyilf 0 -lAf° 

B 2 Tr- 


fr~L( f-Af°;B-|a o -1|f 0 ) J 


■(x;T*/2) 


-(f;B-|a o -l| f 0 ) df 


2 (B-|a -l|f -|Af°|) 

S/N a 2 -2-2- j i_(x;T*/2) 

(B-|ct --iy ) T* 


(Eq. D.3) 
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M 

if the spread of |y| |y'| is small compared with T* . But |y'[ = IYp M ! 
and |y! is just |Yp M I time shifted; thus the spread of |y| |y'| is 
roughly the spread of |Yp M I and will be small compared with T* if the 
spread of lYp^l is small. The magnitude of the fourier transform of 
E[y(x)] is 

l3?[E(y(T))]| = |E[Y(f)]| «E[|Y(f)|] 

« S/N a^(B-|o to -l|f 0 -|Af°|)/(B-|a o -l|f 0 )]sinc( 1 T f T*) . 

(Eq. D.4) 


The sine function is 50% of its peak value at f = ± 0.6 it ; therefore 
E[Y(f)J has 50% of its peak value at f = + 0.6/T* . The slope at 
f - - 0.6/T* is; 


0.43 tt T* S/N a^(B-|cc 0 -l|f 0 -|Af 0 !)/(B-|a o -l|f 0 ) (Eq. 0.5) 

From the above it is clear that if the spread of |Ypj^1 2 is less than T*, 

the echo portion of the matched filter otuput for each of the FM pulses 
is a segment T* in length of band-limited, gaussian, white noise of 

M 

bandwidth B-1ot Q -l[f Q , the bandwidth of Tpp|(f> a 0 ) , while the noise 
portion is a segment Tp in length of gaussian white noise of bandwidth B : 


x(t) « [h 1 (t)j - u(x;T*/2) + n( T )4-u( T ;T F /2)]/ /T f 
x'(T)«[h li (T).n.(T;T*/2) + n'( t)j^(t;T f /2)]//T f 


y(T)w{[h'*(T)h"(T) + h , *(T)n , (x)+n*(T)h l, (i)]j-L/T;T*/2) 
+ n*(T)n'(T)-^(T;T F /2)} /T p . 


(Eq. D.6) 
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The centre frequencies are f g + Af°, f Q - Af° , f g and f respectively 

for h'(x), n(x), h"(x) and n'(x). From the discussion in the main test 
EC|Y(f) - |Y e (f)|] 2 and E{[|Y(-f)| - |Y e (-f)|] [Y(f) - |Y e (f)|]} are 

sought where 


v e ( f ) = ;r[y e ( T )] * 

y e (x) = h'*(x) h 11 (xx,T*/2) 


(Eq. D.7) 


If, as is assumed, ! Y e (f)| is large compared with |Y(f)| - |Y (f)| , 
then |Y(f)| - |Y (f)| is approximately the variation of Y(f) along the 
direction of Y (f) : 


I Y(f) | - | Y (f) | * Re[e'^6Y(f)] , 


(Eq. D.8) 


where <p is the argument of Y g (f) and 6Y(f) = Y(f) - Y g (f) . 
Furthermore, since Y (-f) « Y g *(f) , 


IY(-f) | - | Y (-f) | * Re[e +j<t \5Y(-f)] . 


(Eq. D.9) 


Let <5Y(f) have real and imaginary parts 6Y R (f) and <SYj(f) i then 

E[|<5Y(f)| 2 ] = E{[6Y R (f)] 2 } + E{ [6Y j (f) ] 2 } , 

E{[6Y(f)] 2 } = E{[<5Y R (f)] 2 } - E{[6Yj(f)] 2 } + 2j E{<5Y R (f) 6Yj(f)} , 
E{Re[e' j4, 6Y(f)]} = E{[6Y R (f) cos<i>+<5Yj (f) sin<f>] 2 } 

= 1/2 E[|6Y(f)| 2 ] 


(Eq. D.10) 
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if E{[6Y(f)] 2 } = 0 . Furthermore 

E[<5Y*( -f ) 6Y(f)] = E[6Y R (-f) 6Y R (f)] + E[6Yj(-f) 6Yj(f)] 

- .1 E[6Y R (-f) 6Yj(f) - 6Yj(-f) 6Y R (f)] 

E[6Y(-f) 6Y(f)] = E[5Y R (-f) 6Y R (f)] - E[6Yj(-f) 6Yj(f)] 

- J E[6Y R (-f) <5Yj(f) + 6Yj(-f) <5Y R (f)] 

E{Re[e + ^6Y(-f)] Re[e' J<J 6Y(f) ]} = E{[Y R (-f) cos* - Yj(-f) sin*] 

[Y R (f) cos* + Y(f) sin*]} 

= 0 

(Eq. D.ll) 

if E[6Y*(-f) 6Y(f)] = E[6Y(-f) 6Y(f)] . 

From the properties of the fourier transform, the inverse transforms of 
16Y(f) | 2 , <5Y(f) 2 , 6Y*(-f) 6Y(f) , 6Y(-f) 6Y(f) obey 

jf 1 [ I<SY(f) | 2 ] = 6y*(-T) x 6y(x) , 

97 _1 [6Y(f) 2 ] = <5y(x) x 6y(x) , 

^ _1 [6*Y(-f) 6Y(f)] = 5y*(x) x 6y(x) , 

^ _1 [6Y(-f) 6Y(f)] = 6y(-x) x 6y(x) , 

(Eq. D.12) 

where 6 (x) = y(x) - y (x) ; therefore we need to evaluate the 

expected values of these four convolutions and show that it is zero for 
the last three. The second and last convolution clearly have zero 
expectations since they are sums of terms such as 

E[h'*(x±o) h'*(o)] E[n'(x±o) n'(a)] , 

which are zero. On the other hand, since h'(x), h"(x) and n'(x) are 
uncorrelated: 
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E[y*( -t ) x y(x)] « A(t) + B(r) + C(t) 

E[y*(x) X y(x)] «A-(x) + b'(t) + c'(x) 

where (Eq. ^.13) 

A(x) = E{[n*(-x) h"(-x )j i_(x;T*/ 2)]* ® [n*(x) h"(x)-ro-(x;T*/2)]}/T 2 , 

x 

A'(x) = E{[n(x) h"(z)^-(r;T*/2)]* ®[n*(z) h"(x)j-i-(x;T*/2)]}/T 2 , 

x 

B(x) = E{[h'*(-x)n'(-x)^.(-x;T*/2)]*®[h'*(x)n , (x)-r-v(x;T*/2)]}/T^ , 

x 

B'(x) = E{[h'*(x) n 1 (x) j (x;T*/2)]* ® [h **(x)n 1 (x)^(x;T*/2)]}/Tf , 

x 

C(x) = E{[n*(-x) n'(- x)J~>-(-t;T p /2)]*® [n*(x) n 1 (t)j-l(t;T p/2)]}/T* , 

x 

C (x) = E {[n* (x) n , (x)-r-^(x;Tp/2)]* ®[n*(T) n ' (x).r-u(x;T F /2)]}/T* . 

x 

(Eq. D.14) 

Consider evaluation of C(x) and C‘(x) first, since these are the 
easiest: 


C(x) = /-^(x+o;T F /2) J ~x(a;T F /2) E[n(x+a)*n(a) J E[n'*(-x+a)n 1 (a)] da/T 2 , 


a 00 


— /j-x.(x+o;T f /2)-t-l(o;T f /2) sinc 2 (7rBx)e" 27TjAf ° T da ; 


Tp 


!L sinc 2 (TrBx) e '^ jAf ° T 


(Eq. D.15) 


and 


66 






SACLANTCEN SM-155 


C'(t) = /-n-(x-o; Tj-/2)j-<-{o;T p/2) E[n(x-o)n*(o)] E[n’*(x-o)n 1 (a)]da/Tp , 


0 °° 


— /j-T-(x-a;T F /2)rT.(a;T F /2)sinc 2 [TiB(x-2o)]e" 27TjAf (r ‘ 2a) do , 


T 2 — 


a n -(B-|Af°|) 


2B 2 T 2 


(t,T f ) 


(Eq. D.16) 


since n(x) and n‘(x) are independent and BT F is small compared with 
unity. 

The fourier transforms of C and C 1 are 


j:[c(x)] 


BT r 


(f+Af°;B)(1 -1f+Af°|/B) , 


y[C'(x)] 


On^B-lAf 0 !) 


B ? T r 


sinc(2irfT F ) 


(Eq. D.17) 


In the same manner it may be shown that since BT* is large compared 
with unity 


F [B(x)] « S/N -2- -r-x- [f + (o 0 -l)f 0 /2+Af 0 ;|a 0 -l|f 0 /2 , B-(a Q -l |f Q /2] 
BT F 


rCB'(T)] ^ S/N sinc(27ifT*)^-[(a o -l)f o ^+Af 0 ;|a o -l|f 0 /2,B-|v 1 ' f o /2] 

BT r 


^CA(x)] * S/N -H- Cf-(a 0 -l)f 0 /2+Af 0 ;|a 0 -l|f 0 /2 , B-|a rt -l | f ft /2] 


BT 


0 1 0' 


FCA'(t)] * S/N -H- sinc(2rrfT*)^[-(a 0 -l)f 0 /2 + Af 0 ;|a 0 -l|f 0 /2,B-!a 0 -l|f 0 /2] 
BT F 

(Eq. D.18) 
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where , 

i m ^ f. 

■n.(fs F,, F 2 ) = Jo | f| > F 2 

((F 2 -f)/(F 2 -F 1 ) otherwise . (Eq. D-19) 


With these results E[|6Y(f)J 2 ] and E[6*Y(-f)6Y(f)] are readily 
determined. For f = 0.6/T* they may be further approximated to: 


(S/N) 2 a n " 
E[ | <5Y(f) | 2 ] » - 


B-(a o -Tff 0 -|Af°/2[ 

B-|a -1 1 f 

1 o ' 0 


2 B-|Af°| 1 

— + -- 

S/N B (S/N) 2 


(S/N) 2 a 


B T c 


2 1 
+ ‘ 

S/N (S/N) 2 


(Eq. D.20) 


Since BT* is large compared with unity and sine (1.2n) and 
sine (1.2tt T^/T*) are negligible. The final result is thus Eq. 20, 

since 0.43-rr /T’«2. These results differ from those of Plaisant for 
several reasons: In the first place he assumed Tp = Ax^ , in the 

second place he approximated E{[|Y(f)|-|Y g (f)|] 2 ) by E[|6Y(f)| 2 ] 

instead of by the correct E[|6Y(f)| 2 ]/2. 


D.2 Correlation in time 

The procedure to be followed is the same, except that the roles of time 
and frequency are reversed. Obviously this procedure must yield the same 
approximate variance at half-correlation-peak points since the variation at 
the correlation peak is used; furthermore, the covariance between the 
half-correlation-peak points is again approximately zero. However, in the 
expected correlation function, sine (uf T*) must be replaced by 

sine [tit (B-1 ot Q - 11f Q -1Af°| ], since B-1 ct o ~l |f Q is the analog of T* ; 

therefore, in the slope, T* must be replaced by B - i a 0 ~l!f 0 ancl ^ becomes 


0.43 tt S/N a 2 (B-|a o -l|f 0 -|Af 0 |) 2 /(B-| % -l|f 0 ) 

The final result is thus Eq. 0.21. 


(Eq. 0.21) 
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D.3 Phase of the Peak 

To measure the phase at the maximum magnitude of double correlation 
it is necessary only to take the arctangent of the ratio of the imaginary 
part to the real part of the double correlation at this point. The error 
in measuring phase is just the variation in correlation perpendicular to 
the expected correlation divided by the magnitude of the expected 
correlation. But the variance of the correlation in this direction is, 
following the reasoning of D.l, just one-half the total variance of the 
correlation E[|6Y(f)| 2 ] , which is given at the end of Sect. D.l. If 
the square root of half this quantity is divided by the magnitude of 
the expected maximum correlation, which from the beginning of Sect. D.l 
is S/N o n 2 (B-|a Q -l|f 0 -|Af°|)/(B-|a o -ljf Q ) , the result is Eq. 44 

of the main text. 
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